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Abstract 

in 

Cn In this paper we present an overview of pattern formation analysis for an analogue of the Swift- 

^.^ Hohenberg equation posed on the real hyperbolic space of dimension two, which we identify with the 

^S^ Poincare disc D. Different types of patterns are considered: spatially periodic stationary solutions, 

"T"^ radial solutions and travehng waves, however there are significant differences in the results with the 

,S^ Euclidean case. We apply equivariant bifurcation theory to the study of spatially periodic solutions on 

5^ a given lattice of D also called H-planforms in reference with the "planforms" introduced for pattern 

C formation in Euclidean space. We consider in details the case of the regular octagonal lattice and give 

' ' a complete descriptions of all H-planforms bifurcating in this case. For radial solutions (in geodesic 

^vj polar coordinates), we present a result of existence for stationary localized radial solutions, which we 

^ have adapted from techniques on the Euclidean plane. Finally, we show that unlike the Euclidean case, 

I^Tj the Swift-Hohenberg equation in the hyperbolic plane undergoes a Hopf bifurcation to traveling waves 

A—^ which are invariant along horocycles of D and periodic in the "transverse" direction. We highlight our 

\^ theoretical results with a selection of numerical simulations. 

f~^ Keywords: Swift-Hohenberg equation; Pattern formation; Poincare disk; Equivariant bifurcation; Traveling 

f^ wave. 
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^ 1 Introduction 

The goal of this paper is to lay the foundations of pattern formation theory for PDE's which are defined in 
an hyperbolic space and more specifically in the hyperboUc plane, which we identify from now on with the 
Poincare disc D. The origin of our interest in this question comes from a model for the visual perception of 
textures by the cortex of mammals, in which it is assumed that textures are represented by the "structure 
tensor" , a 2 x 2 symmetric, positive definite matrix, and the neurons in the visual cortex region named 
VI are sensitive to the values of the structure tensor. It can be shown that there is no loss of generality 
to assume that the structure tensor has determinant 1, and it turns out that the space of such structure 
tensors is isomorphic to the hyperboUc plane. Spontaneous activity in the cortex then leads to a bifurcation 
problem with pattern formation. Several papers have been published on this topic [17, 18, 29-31]. The model 
equations are integro-differential. However under certain conditions they can be converted to PDE's which 
look like the Swift-Hohenberg equation in E^, except that the Laplacian is replaced by the Laplace-Beltrami 



operator Ad in D [28]. We shall therefore concentrate on this equation now defined in D, which we take of 
the general form 

ut = -{a'^ + Anfu + Xu+Af{u), z G D, (1.1) 

where 

7V(u) = lyu^ - riu^ 

and A, a, v and r\ are real coefficients. There is another reason to consider (1.1). The Swift-Hohenberg 
equation was initially derived to model pattern formation in the Benard problem, that is the onset of 
convection for a fluid flow set between two horizontal and extended plates with a positive temperature 
gradient across the shell [63] . This equation turns out to be the " simplest" PDE that contains the ingredients 
for steady state pattern formation from an homogeneous state: symmetries (Euclidean invariance), existence 
of a critical non zero wave number, nonlinearity. Most basic phenomena that are associated with pattern 
formation in physical and chemical/biological systems are captured by this equation. 

We may similarly consider the Swift-Hohenberg equation in hyperbolic space (or plane) as a model equation 
for pattern and wave formation in problems related to quantum chaos [4, 7, 23, 61] and cosmological theories 
[22, 40, 50]. Note that more recently, some aspects of dispersive and concentration phenomena have been 
studied for evolution equations such as the Nonlinear Schrodinger equation and the wave equation posed on 
hyperbolic space [2, 3, 8]. The methods presented here should be transposable to these types of equations. 

Our aim is to "translate" to the hyperbolic case the methods which have been successful in the Euclidean 
case. As we said above, these methods are not restricted to the Swift-Hohenberg equation and can be applied 
to many other systems. However this equation is sufficiently simple and in a sense, "universal", to serve as 
a guideline for our study. 

In order to facilitate the reading of the paper we first recall below the main ideas, which lie the analysis of 
pattern formation in Euclidean space. 

The mathematical modeling of pattern formation in morphogenesis originates in the seminal work of Alan 
Turing on "The Chemical Basis of Morphogenesis" [65] (1952). The model is based on systems of reaction- 
diffusion equations. A theory of pattern formation in physical systems was developed at about the same 
time by Lev Landau who derived and analyzed with V. L. Ginzburg an equation for phase transition in 
supraconductivity, which turned out to be adequate in many other problems involving phase transitions [48] . 
Another root of pattern formation theory can be found in the work of G.I. Taylor on stability of Couette flow 
(primary flow of a fluid filling a cylindrical shell with rotating inner cylinder) [64] (1923). There the flow is 
governed by Navier-Stokes equations, as is the Benard problem of onset of convection. In the framework of 
pattern formation all these problems share similar properties with the Swift -Hohenberg equation, which we 
may formally write 

^=F(.,A) (1.2) 

where u = u{x,t) e M is a scalar field, x E M" (in general ?i = 1, 2 or 3), and i^ is a smooth operator 
defined in some Banach space. The coefficients a, v and t] are fixed. This equation is invariant under 
isometric transformations of U: let us define Tgu{x,t) — u{g~^ ■ x,t) where g is any isometry in M". Then 
F{TgU, A) = TgF{u, A) for all g € E(n, M) (EucHdean group), x € M" and t G K. The linear stabihty analysis 
of the trivial, homogeneous, state u = is done by solving the linearized equation 

ut = —{o? -t- Ak")^u + \u 

for disturbances of the form u(a;, t) — e'^*(ylcos (k • x) + S sin (k • a;)) as A is varied. The vectors k G M" are 
called wave vectors. This leads to a dispersion relation 

a = -i-\\k\\^ + a^f + X. 

When A < all such modes are exponentially damped to as i — > +oo^. At the critical value A = the 
eigenvalues a fill the half-line (— oo, 0]. For the critical eigenvalue the "neutral" modes are A cos (k ■ x) + 

^It can be shown that is globally attracting against periodic perturbations for the Swift-Hohenberg equation thanks to its 
variational structure. 



B sin (k • x) where k is any vector with ||k|| = a. This value a is the critical wave number. When A > the 
critical modes become unstable and we may expect the bifurcation of solutions which are not homogeneous 
but instead have a wavy structure in space. However a continuum of Fourier modes become simultaneously 
unstable. Moreover, to the critical wave number there corresponds a full circle of critical wave vectors with 
that length. These facts forbid us from directly applying standard reduction techniques such as Lyapunov- 
Schmidt decomposition [15] and center manifold theorem [36], to the bifurcation analysis of this instability. 
Indeed, a crucial hypothesis in these methods is that the neutral (or center) part of the spectrum of the 
critical linearized operator can be separated from the rest of the spectrum and moreover consists of a finite 
number of eigenvalues with finite multiplicity. We therefore need to be more specific in the type of patterns 
we are looking for. 

There are basically two main cases which have been thoroughly investigated and which lead to elegant as well 
as relevant solutions from the observational point of view: (i) periodic patterns and (ii) radially symmetric 
(localized) patterns. Other types of patterns have also been studied by various authors, such as quasi periodic 
patterns (see [43] for recent advances in this case), spirals, defects in periodic patterns and others. All of 
these are important from physical point of view. However we concentrate in this paper on the fundamental 
patterns (i) and (ii) . 

Let us briefly recall the basic ideas which lie behind the two main approaches. 

(i) Periodic patterns. In the 60's and early 70's several authors, among whom Klaus Kirchgassner 
played a prominent role, established a rigorous nonlinear theory for the bifurcation of cellular solutions 
in the Couette- Taylor and Benard problem, see [47] for a general exposition. The idea was to restrict 
the problem to spatially periodic solutions in order to obtain a finite dimensional bifurcation equation 
by applying Lyapunov-Schmidt decomposition for Navier-Stokes equations. Later these ideas have 
been generalized and geometrically formalized (see [15] for a bibliography). On the circle of critical 
wave numbers we can select any two non colinear vectors ki and k2, and span a periodic lattice 
C* = {n-iki + n2k2 ,ni,n2 € Z}. Any harmonic function exp (ik • x) with k € £*, is biperiodic in the 
plane with respect to translations piei +P2G2, where e^ e M^ are such that e.^ • kj — 6ij and pi, p2 are 
integers. The set of these translations forms a lattice group C (£* is the dual lattice of C). Suppose now 
we look for the bifurcation of solutions which are invariant under the action of C. These solutions are 
therefore spatially periodic. This implies that the equation can be restricted to the class of £ periodic 
functions and the symmetry group of translations M^ has now to be taken "modulo translations" in £, 
in other words the translation group of symmetries of the equation becomes M^/Z^ ~ T^ (the 2-torus). 
Moreover within this class, the spectrum of the linearized operator is discrete with finite multiplicity 
eigenvalues. In particular the critical eigenmodes, which corresponds to those harmonic functions in 
C with ||k|| = a, are in finite number. This allows to apply the center manifold reduction theorem 
and get a bifurcation equation in M.'^ where d is the number of wave vectors in C* with critical length. 
There are three possible cases: 

1. If ki_Lk2, £ is a square lattice. Then d = 4 (critical wave numbers ±ki, ±k2). C is invariant 
under rotation of angle tt/2 as well as under reflections through the axis along ei for example. 
This spans the symmetry group of the square D4. The full symmetry group of the equation in 
this class of functions is the semi-direct product D4 k T^ . 

2. If ki and k2 belong to the vertices of a regular hexagon, £ is a hexagonal lattice. Then d = 6. C 
is invariant under the symmetries of the hexagon, which form the 12 element group Dg. The full 
symmetry group of the equation in this class of functions is the semi-direct product Dg k T^ . 

3. In all other cases the lattice is rhombic (rectangular). This however implies that there exist in 
C* vectors with length smaller than kc. Such a vector can be chosen as a basis element for the 
lattice. However perturbations with the corresponding wave function are damped and should 
not be observable in general (certain conditions may allow for these modes to coexist as critical 
modes, we do not enter in these cases here). 



Finally, the symmetries of the resulting, finite dimensional system can be exploited to classify the 
bifurcated branches by their isotropy (residual symmetry) and to determine their stability. Details 
can be found in a number of references, see [35], [15] for the general context of equivariant bifurcation 
theory, and [39] for a thorough exposition of the mathematical theory of pattern formation. 

It turns out that in hyperbolic geometry the same procedure applies but leads to very different results. 
The reason is that there are infinitely many types of lattices in D. Given an integer n > 3 one can 
always build an n-gon in D which tiles the entire hyperbolic plane, while in contrast there are only 4 
types of fundamental tiles in M^ (oblique, rhombic, square and hexagonal). Moreover if F is a lattice 
group in D (we give a precise definition in Section3), then the compact surface D/F is a torus, the genus 
of which can take any value p > 2 (depending on F). This torus admits a finite group of automorphisms 
(again in contrast to the Euclidean case where the symmetry group of the torus has the form D^. x T^ 
with k ^ 2 (rhombic), 4 (square) or 6 (hexagonal). This will be exposed in Section 3. We have applied 
the method to the simplest case of a regular octagonal lattice. All possible bifurcations are described 
and it is shown that non trivial dynamics can also bifurcate generically in one case. 

(ii) Localized states and radially symmetric patterns. Experiments in pattern formation can also 
produce localized solutions: fronts, bumps, spirals, to name a few. The method exposed above for 
analyzing the bifurcation of periodic patterns does not work and another approach must be undertaken. 
The fundamental idea was introduced in the early 80's by Klaus Kirchgassner [46] in the context of 
solitary water waves (Korteweg-de Vries equation) and further developed by G. looss and others (see 
[32]). Suppose for a moment that the domain is one dimensional: a: G M and one is looking for steady 
states: F{u,X) = 0. This equation is an ODE in the space variable a;: = —{a^+d'^)^u+Xu+i'u^—r]u^. 
Bounded solutions of this ODE will produce admissible patterns. This equation can be written as a 4 
dimensional evolution equation 
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where U — u,V — dxU, W — d^u and Z = d^u. The linear part has a pair of purely imaginary, double 
eigenvalues ±ja at A = 0. Another essential property of this system is "time" reversibility: changing 
X to — cc, V to —V and Z to —Z does not change the equations. This bifurcation problem corresponds 
to a reversible Hopf bifurcation with 1 : 1 resonance, which has been studied in details for example 
in [41, 42]. Applying reversible normal form theory (see [36, 41, 42]) one recovers the bifurcation of x 
periodic patterns, but also one can show depending on the values of coefficients i/ and j] the bifurcation 
of homoclinic orbits to the trivial state or to periodic orbits, which correspond to "bumps" decreasing 
to or approaching asymptotically periodic patterns at a; — ^ ±oo. 

Extending this approach to multidimensional problems is technically involved. A rigorous theory was 
handled in the 90's by Scheel [59, 60] for the bifurcation of radially symmetric solutions of reaction- 
diffusion equations. His idea was again to consider the system as an evolution problem, but now 
with respect to the radial variable r. This leads to a singular problem which can be tackled by 
applying successive changes of variables and a suitable center manifold theory. A detailed analytical 
and numerical investigation has also been presented in [6, 11-13, 51-54] for the Swift-Hohenberg 
equation. We shall detail the procedure in the context of hyperbolic geometry in Section 4 and show 
that in some respect, the situation is simpler in hyperbolic than in Euclidean plane. 

A final remark in this introduction is that pattern formation has been intensively studied in Euclidean and 
spherical geometry, see [16, 24, 25], and the present study completes the landscape by adding hyperbolic 
geometry to the classification. We believe this gives an intrinsic interest to this study. 

In the next section we expose basic notions of the geometry and harmonic analysis in the Poincare disc, in 
particular the Fourier Helgason transform which is the equivalent to Fourier transform in M^. 



2 Geometry and harmonic analysis in the hyperbohc plane 

2.1 Isometries of the Poincare disk 

The hyperbolic plane is a Riemannian space of dimension 2 with constant curvature —1. It admits several 
representations among which the most suitable for our purpose is the "Poincare disc" D = {z G C | |z| < 1} 
equipped with the distance 



do{z, z') ^ 2atanh ( ' ^ ' j . (2.1) 

Geodesies in D are carried by circles which intersect the boundary 9D orthogonally. The corresponding 
measure element is given by 



dm(^) = 7^^^^- (2-2) 



Adzdz 

We now describe the isometries of D, i.e the transformations that preserve the distance do. We refer to 
the classical textbook in hyperbolic geometry for details, e.g, [45]. The direct isometries (preserving the 
orientation) in D are the elements of the special unitary group, noted SU(1, 1), of 2 x 2 Hermitian matrices 
with determinant equal to 1. This is a 3-parameter real simple Lie group. Given: 

7 = ( - ^1 such that lap - |/3|2 = 1, (2.3) 

\ P a ) 

an element of SU(1, 1), the corresponding isometry 7 in D is defined by: 

7-^=f^, ^eD. (2.4) 

pz -\- a 

Orientation reversing isometries of D are obtained by composing any transformation (2.4) with the reflexion 
K : z — >■ z. The full symmetry group of the Poincare disc is therefore: 

U(l, 1) = SU(1, 1) U K • SU(1, 1). 

Let us now describe the different kinds of direct isometries acting in D. We first define the following one 
parameter subgroups of SU(1, 1): 

^a^ JcoMr/2) Binh(r/2)\ 

^ \^ sinh(T/2) cosh(T/2) / ^ 

n"^ {n,= I : ' ^" ^ ^" . ),seM}. 




Note that r^ • z = e"^z for z G D and also a^- • = tanh(T/2). The following decomposition holds (see [44]). 
Theorem 2.1 (Iwasawa decomposition). 

SU{1, 1) = KAN. 

The group K is the orthogonal group SO (2) which fixes the center O of D. Its orbits are concentric circles. 
The orbits of A converge to the same limit points of the unit circle 9D, fej-i = ±1 when r — > ±00. The 
elements of A are sometimes called boosts in the theoretical Physics literature [7] . They are circular arcs in 
D going through the points 61 and 6_i. The orbits of TV are the circles inside D and tangent to the unit 
circle at 61. These circles are called horocycles with base point 61. TV is called the horocyclic group. These 
orbits are shown in Figure 1. 

The Iwasawa decomposition allows us to decompose any isometry of D as the product of at most three 
elements in the groups, K, A and N. Then, it is possible to express each point z e D as an image of the 
origin O by some combination of elements in K, A and N. There are two important systems of coordinates: 




bi b-i 




Figure 1: The orbits in the Poincare disk of the three groups K, A and N. 

1. the geodesic polar coordinates: r^a^ • O = tanh(T/2)e*'^ with r = do{z,0), 

2. the horocyclic coordinates: z — Usar • O € D, where rig are the transformations associated with the 
group N (s S M) and a^. the transformations associated with the subgroup A (r G M). 



2.2 Periodic lattices of the Poincare disk 

A Fuchsian group is a discrete subgroup F of SU(1, 1). We are going to be concerned with fundamental 
regions of Fuchsian groups. 

Definition 2.1. To any Fuchsian group we can associate a fundamental region which is the closure, noted 

o 

Fr, of an open set F-p C D with the following properties: 

(i) if-f^ Ide F, then 7 • Fr n i^r = 0; 
(it) U 7 • ^r = D. 

The familly {^ ■ F-p \ j £ T} is a tesselation o/D. 

Fundamental regions may be unnecessarily complicated, in particular they may not be connected. An 
alternative definition is that of a Dirichlet region of a Fuchsian group. 

Definition 2.2. Let T he a Fuchsian group and z G D &e not fixed by any element ofT\ Id. We define the 
Dirichlet region for F centered at z to he the set: 

D,(F) = {z'eD I dn{z\z) < do{z n ■ z)^! eT}. 

From [45], we have the following theorem. 

Theorem 2.2. If z (zV is not fixed by any element ofT\Id, then Dz{T) is a connected fundamental region 
forV. 

Let F be a Fuchsian group acting on D with /i(D/F) < 00, and F^ be a fundamental region for this action. 
We write tt : D — > D/F the natural projection and the points of D/F are identified with the F-orbits. The 
restriction of tt to Fp identifies the congruent points of Fr that necessarily belong to its boundary 9Fr, 
and makes D/F into an oriented surface. Its topological type is determined by the number of cusps and by 
its genus: the number of handles if we view the surface as a sphere with handles. By choosing Fr to be 
Dirichlet region, we can find the topological type of D/F (in this case D/F is homcomorphic to Fp/F, see 
[45]). Furthermore, if finite, the area of a fundamental region (with nice boundary) is a numerical invariant of 
the group F. Since the area of the quotient space D/F is induced by the hyperbolic area on D, the hyperbolic 
area of D/F, denoted by /^(D/F), is well defined and equal to /x(Fr) for any fundamental region Fp. If F has 
a compact Dirichlet region Fr, then Fp has finitely many sides and the quotient space D/F is compact. If, 
in addition, F acts on D without fixed points, D/F is a compact Riemann surface and its fundamental group 
is isomorphic to F. 



Definition 2.3. A Fuchsian group is called cocompact if the quotient space D/F is compact. 

When a Fuchsian group is cocompact, then it contains no parabolic elements and its area is finite [45]. 
Furthermore a fundamental region can always be built as a polygon. The following definition is just a 
translation to the hyperbolic plane of the definition of an Euclidean lattice. 

Definition 2.4. A lattice group o/D is a cocompact Fuchsian group which contains no elliptic element. 

The action of a lattice group has no fixed point, therefore the quotient surface D/F is a (compact) manifold 
and it is in fact a Riemann surface. A remarkable theorem states that any compact Riemann surface is 
isomorphic to a lattice fundamental domain of D if and only if it has genus g > 2 [45]. The case .g = 1 
corresponds to lattices in the Euclidean plane (in this case there are three kinds of fundamental domains: 
rectangles, squares and hexagons). The simplest lattice in D, with genus 2, is generated by an octagon and 
will be studied in detail in Section 3. 

Given a lattice, we may ask what is the symmetry group of the fundamental domain Fy, identified with the 
quotient surface D/F. Indeed, this information will play a fundamental role in the subsequent bifurcation 
analysis. In the case of Euclidean lattice, we recall that the quotient M^/F is a torus T^ (genus one surface), 
and the group of automorphisms is 'H ix T where H is the holohedry of the lattice: H = D2, D4 or Dg for the 
rectangle, square and hexagonal lattices respectively. In the hyperbolic case the group of automorphisms of 
the surface is finite. In order to build this group we need first to introduce some additional definitions. 

Tilings of the hyperbolic plane can be generated by reflections through the edges of a triangle T with vertices 
P, Q, R and angles n/i, tt/tu and ir/n respectively, where £, m, n are integers such that l/£+l/m + l/n < 1 
[45]. 




Figure 2: The triangle T. The values of I, m, and n are I — 8, m — 2 and n = 3. 

Remember that reflections are orientation-reversing isometrics. We note k, k' and k" the reflections through 
the edges PQ, QR and RP respectively (Figure 2). The group generated by these reflections contains 
an index 2 Fuchsian subgroup A called a triangle group, which always contains elliptic elements because 
the product of the reflections through two adjacent edges of a polygon is elliptic with fixed point at the 
corresponding vertex. One easily shows that A is generated by the rotations of angles 27t/1, 27r/m and 2TT/n 
around the vertices P, Q, R respectively. A fundamental domain of A is the "quadrangle" Fa = t U kt 
[45]. Note that Fa — D/A is a sphere (genus surface) obtained by identifying the three edges of T. The 
subgroup of hyperbolic translations in A is a lattice group F, normal in A, whose fundamental domain is 
flUed with copies of the basic tile T. The group of orientation-preserving automorphisms of Fp ~ D/F is 
therefore G = A/F. From the algebraic point of view, G is generated by three elements a, b, c satisfy the 
relations a^ — fe™ — c" — I and a ■ b ■ c ^ 1. We say that G is an {l,rn,n) group. Taking account of 
orientation-reversing isometries, the full symmetry group oi Fr is G — G U nG — G x TL^i^). This is also a 
tiling group of Fp with tile T: the orbit Q ■ T flUs Fp and its elements can only intersect at their edges. 

Given a lattice, how to determine the groups G and Ql The following theorem gives conditions for this, see 

[37]. 

Theorem 2.3. An (l, m, n) group G is the tiling rotation group of a compact Riemann surface of genus g if 
and only if its order satisfies the Riemann- Hurwitz relation 

\C\- ^9^^ 



Tables of triangle groups for surfaces of genus up to 13 can be found in [10]. 

2.3 Laplace-Beltrami operator 

2.3.1 Definition 

The Laplace-Beltrami operator in D is defined by the expression in cartesian coordinates z ^ zi + iz2 G D: 

It has the fundamental property of being equivariant under isometric transformations. More precisely, let 
us define the following transformation on functions in D. We set 

Tgfiz) = /(5-1 • z) for all g G U(l, 1), z e D. 

Then AoTgf = TgAof. The set of transformations Tg defines a representation of the group U(l, 1) in the 
space of functions in D. 

2.3.2 General eigenfunctions of the Laplace-Beltrami operator 

Let 6 be a point on the circle 9D. For z e D, we define the "inner product" (z, b) as the algebraic distance 
to the origin of the (unique) horocycle based at b and passing through z. This distance is defined as the 
hyperbolic signed length of the segment O^ where ^ is the intersection point of the horocycle and the line 
(geodesic) Ob. This is illustrated in Figure 3 in the case b = bi = 1. Note that {z,b) does not depend on 
the position of z on the horocycle. In other words, {z, b) is invariant under the action of the one-parameter 
group N. Each point z of D can be written in horocyclic coordinates z = n^ ■ {a-^ ■ O). Then (z,fei) = t. 
Note that r is negative if O is inside the horocycle and positive otherwise. 




Figure 3: The construction of the "signed inner product" (z, b) with b — 1. 

In analogy to the Euclidean plane waves, we define the "hyperbolic plane waves" as the function 

ep,fc(z) = eC"'+5)<-'^>, peC. (2.6) 



Lemma 2.1. For the Laplace- Beltrami operator defined in equation (2.5), we have 

Ad ep,b{z) = - (p' + ^ j ep,fc(z), V {p, b,z) eCxdDxD. (2.7) 

Proof. See [38]. D 

These elementary eigenfunctions allow the construction of general eigenfunctions of Ad. Let A{dD) denote 
the space of analytic functions on the boundary dD of the Poincare disk, considered as an analytical manifold. 
Let U be an open annulus containing 9D, 'HiU) the space of holomorphic functions on U equipped with the 
topology of uniform convergence on compact subsets. We identify A{dD) with the union [JijHiU) and give 
it the limit topology. The element of the dual space A'{dfii) are called analytic Junctionals or hyperf unctions. 
Since elements oi A'{dli) generalize measures, it is convenient to write 



m(/) = /" f{b)dxn{b), / e yl(aD) and m e yl' 

Jdn 



From Helgason's theory [38] we have the following theorem. 

Theorem 2.4. The eigenfunctions of the Laplace- Beltrami operator on D are precisely the functions 



*(z) = / ep4z)dmp{b), (2.8) 

JdB 
where p £ C, nip £ A'{dD) and the eigenvalue is — (p^ + j) . 

Note that real eigenvalues — (p^ + j) of Ad correspond to taking p real or p £ iR. The latter case is 
irrelevant for the following study as it corresponds to exponentially diverging eigenfunctions. Therefore the 
real spectrum of Lb is continuous and is bounded from above by — | . 

2.3.3 Periodic eigenfunctions of the Laplace-Beltrami operator 

In the following we will look for solutions of bifurcation problems in D, which are invariant under the action 
of a lattice group: 7 • u{z) — u{j^^ ■ z) — u{z) for 7 e F. This reduces to look at the problem restricted to 
a fundamental region with suitable boundary conditions imposed by the F-periodicity, or, equivalently, to 
looking for the solutions of the problem projected onto the orbit space D/F (which inherits a Riemannian 
structure from D). Because the fundamental region is compact, it follows from general spectral theory that 
Ad is self-adjoint, non negative and has compact resolvent in L^(D/F) [14]. Hence its spectrum consists of 
real positive and isolated eigenvalues of finite multiplicity. 

Coming back to Theorem 2.4, we observe that those eigenvalues A of Ad which correspond to F-invariant 
eigenfunctions, must have p S M or p e iR. The case p real corresponds to the Euclidean situation of 
planar waves with a given wave number, the role of which is played by p in D. In this case the eigenvalues 
of Ad satisfy | < A. On the other hand there is no Euclidean equivalent of the case p £ zM, for which 
the eigenvalues < A < | are in finite number. It turns out that such "exceptional" eigenvalues do not 
occur for "simple" groups such as the octagonal group to be considered in more details in the Section 3. 
This follows from formulas which give lower bounds for these eigenvalues. Let us give two examples of such 
estimates (derived by Buser [14], see also [44]): (i) if g is the genus of the surface D/F, there are at most 
3g — 2 exceptional eigenvalues; (ii) if d is the diameter of the fundamental region, then the smallest (non 
zero) eigenvalue is bounded from below by (An sinh |) 

Suppose now that the eigenfunction in Theorem 2.4 is F-periodic. Then the distribution nip satisfies the 
following equivariance relation [58]. Let 7(0) denote the image oi 9 £ 9D under the action of 7 G F. Then 

mp(7.0) = |7'(0)|^+"'mp(0). (2.9) 

Remark 2.1. As observed by [62], this condition is not compatible with m^ being a "nice" function. In fact, 
not only does there exist no explicit formula for these eigenfunctions, but their approximate computation is 
itself an uneasy task. We shall come back to this point in the next chapter. 



2.4 The Helgason-Fourier transform 

Based on the elementary eigenf unctions 2.6, Helgason built a Fourier transform theory for the Poincare disc, 
see [38] which we recall now. 

Definition 2.5. /// is a complex-valued function on D, its Helgason-Fourier transform is defined by 

/(P, b)^ f /(z)e(-'''+^)<-^''> dmiz) (2.10) 

■/D 

for all p £ C, b £ 9D for which this integral exists. 

If we denote T>{D), the set of differentiable functions of compact support then the following inversion theorem 
holds. 

Theorem 2.5. If f £ V{D), then 

/W = ^ / / fip,b)ei'P+^^y'-''^ptanh{7rp)dpdb (2.11) 

where db is the circular measure on dD normalized by J^p db = 1. 

2.5 Existence, uniqueness and regularity of the solutions 

If we think of u, solution of the Swift-Hohenberg equation (1.1), as a function evolving in space and time 
(for example an averaged potential membrane [18, 31]) it is natural to impose that u should be uniformly 
bounded, i.e., a function in L°°{D). Indeed, the space of bounded functions contains all types of solutions of 
the Swift-Hohenberg: fronts (solutions connecting two homogeneous states) , periodic solutions and localized 
solutions. Note that _LP(D)-spaces with p < oo are only relevant for localized solutions. In addition, this 
natural property should be preserved under the time evolution. In other words, the problem is to know if 
the Swift-Hohenberg equation (1.1) defines a regular semi-flow in the space L°°{D) which is global in time. 
Note first, that the Cauchy problem in this space is not difficult to handle and it is easy to prove that there 
exists an unique solution of (1.1) for a small time interval which depends upon the L°°(D) norm of the initial 
data. In the Euclidean case, this global time problem has received much interests in the past decades and 
can be solved using local energy estimates [20, 21]. The proof of an analog global time existence theorem 
for the hyperbolic plane is far beyond the scope of this paper and instead of working in a functional space 
that suits for all types of solutions, we prefer to use specific functional space for each type of solutions. 

The types of solutions that we are interested in are periodic and localized solutions, which can be defined 
on well chosen Sobolev spaces where the existence of global solutions is straightforward. 

We also recall that the Swift-Hohenberg equation (1.1) posed on D is a gradient system, 

ut = -V£(u), 
in il^(D), where the energy functional £ is given by 

[{a + Ao)u] Xu^ vv? rju'^ 



£iu) 



dm(z), z e D, (2.12) 



and the gradient V£{u) — jf (u) of £ with respect to u is computed in L^(D). In particular, £ decreases 
strictly in time along solution of (1.1) unless the solution is stationary. While we cannot evaluate the energy 
functional along periodic patterns as they are not localized, whence the integral (2.12) may not exist, we 
may, however, define a local energy by integrating over one spatial period of the underlying periodic pattern. 
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3 Bifurcation of F-periodic patterns 



3.1 Reduction of the problem 

3.1.1 Linear stability analysis in the class of F-periodic perturbations 

Let us study the stabihty of the trivial state m = of equation 1.1. If we look for perturbations in the form 
of hyperbolic plane waves e'^*ep^f,(z), solving the linearized Swift-Hohenberg equation comes back to solving 
a "dispersion relation" 

(3.1) 



1. 



Writing z in horocyclic coordinates {s^t) (see 2.3) we may want to look for solutions which are independent 
of the coordinate s along the horocycles based at the point h on 9D and which are periodic in t. This would 
imply p = a + i/2 with a € K. Indeed ep.b(z) = e^'""'"^/^^'^. However in general cr is a complex eigenvalue in 
this case, which we shall treat in Section 5. 

We now look for solutions within the restricted class of functions which are invariant under the action 
of a lattice group T and which are L in the fundamental domain of periodicity. This comes back to 
looking for Equations (1.1) projected onto the compact Riemann surface D/F. As noticed in 2.3.3, the 
spectrum of the Laplace-Beltrami operator in D/F is made of isolated, real and positive eigenvalues with 
finite multiplicities, which we label by increasing order and which we write p^ + |, A: e N. Moreover the 
corresponding eigenfunctions are expressed as in (2.8) with the measure rup^ satisfying the equivariance 
property 2.9. We assume that all eigenvalues are greater than |;, which implies pk € K. This is known to be 
true for example when F is the regular octagonal lattice, a case that we shall analyze in detail in the next 
subsections. 

It follows from the dispersion relation (3.1) restricted to F-invariant perturbations that 

1. When A < all F periodic perturbations are damped and is a stable state. 

2. The neutral stability curve consists in a discrete set of points (p^, A^) on the curve A = [a^ ^ P^ ^ \Y ■ 
There exists a minimum of A^ at some value pk^ as we see in Figure 4. Note that by adjusting the 
value of a we can set A^^ — 0, which we shall assume in the rest of this section. 



Pk 



Figure 4: The neutral stability curve (schematic). 



3. In general the critical parameter value Afe is associated with a unique wave number pk^ ■ In this case 
the critical eigenspace is an irreducible representation of the group Q of symmetries of D/F. 
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3.1.2 Center manifold reduction 

From these points it follows that we can apply a center manifold reduction to equation (1.1). We quickly recall 
the procedure. Let V be the critical eigenspace and x & V. We write V'^ for the orthogonal complement of 
V in L^(D/r). The center manifold theorem asserts the following [15, 36]: 

(i) There exists a neighborhood U of (0, 0) in F x M and a C^ map * : V x R ^ V-^ {p > 2) such 
that the graph of ^ in _L^(D/r) is a locally flow invariant, attracting C^ manifold of equation (1.1) for 
(a;,A) &U. 

(ii) The map ^ is equivariant under the action of the group Q: for any g e Q, '^'{g ■ x,X) — g ■ ^(x, A) for 
aU {x, A) e U. 

Therefore the bifurcation analysis can be reduced to the projection of equation (1.1) onto V with u = 
X + '^{x,X). We write this equation 

f.F(.,A, ,3.2) 

and by point (ii) above F{g ■ x,X) = g ■ F{x, A) for all g S ^ and {x, X) G V x M.. From the assumptions we 
have F{0, 0) = and D^F^O, 0) = 0. The Taylor expansion of F can be computed by a recursive method 
[36]. However V can only be known approximately because the eigenf unctions of Ad do not have explicit 
expressions and are in general difficult to compute (see below). 

We can now apply to (3.2) the methods of equivariant bifurcation theory [15] in order to determine the 
bifurcation diagram. 

3.1.3 Equivariant branching lemma 

When a differential equation like (3.2) is invariant under the action of a symmetry group, this imposes 
geometrical constraints which can be exploited to determine the bifurcation diagram or more generally, 
to analyze the dynamics. Here we recall some basic facts and methods which will be applied in the next 
subsection. 

Definition 3.1. Given a representation (a linear action) of the group Q in a space V , we say that it is 
absolutely irreducible if the only endoniorphisms of V which commute with this action are scalar multiples 
of the identity. 

One can easily check that an absolutely irreducible representation is irreducible, meaning that the only 
subspaces of V which are invaraint under the representation are {0} and V, but the converse may not be 
true. 

We now assume that Q acts absolutely irreducibly in V. This has two consequences: 

(i) F{Q, A) = for all A. Indeed F{0, A) is fixed by all g E Q and the only point which has this property 
is 0. 

(ii) DxF{0,X) — c{X)I (identity operator). 

In the case of equation (1.1), (ii) is automatically satisfied. Indeed the linear part is — (a^ -\- iS.-of'u + Xu and 
we have seen that V is the kernel of —{a^ + A^)^ in iv^(D/r). Let P be the orthogonal projection on V . 
Setting X = Pu, we see that DxF{0, X)x — P{—{a^ + Ad)^^ + Xx) = Xx. 

Definition 3.2. An isotropy group for the action of Q inV is the largest subgroup which fixes a point in V . 
For example Q itself is an isotropy group (it fixes Oj. 
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Different points in V may have the same isotropy group. Let H be an isotropy group. Since the action 
is hnear, the set Vh = {x Cz V \ H ■ x = x} is a subspace of V caUed the fixed point subspace of H . Vh 
contains points with higher isotropy (at least it contains 0) and this induces a stratification of the fixed point 
subspace. The subset of points which have exactly H as isotropy subgroup is open in Vh and is called the 
principal stratum of Vh- The next proposition is straightforward but of high consequences. 

Proposition 3.1. Fixed point subspaces are invariant under equivariant maps in V . 

We can now state the main result of this subsection (equivariant branching lemma for the Swift-Hohenberg 
equation). 

Theorem 3.1. Under the above hypotheses, suppose that H is an isotropy subgroup such that dim Vjf = 1. 
Then a branch ofT-periodic steady states of (1-1) bifurcates from at X — in Vh- Let N{H) be the 
normalizer of H in Q. If N(H)/H ~ Z2, then the bifurcation is a pitchfork in Vh - 

Proof. By the previous proposition, the bifurcation equation (3.2) restricts to the invariant axis Vh- We 
write this scalar equation x = f{x, A) where by assumptions, /(O, A) = and /'(O, A) = A. It follows that we 
can rewrite this equation 

x = x[\ + h[x,\)] (3.3) 

where h{x. A) — 0{x^). Having eliminated the solution a; == we obtain the bifurcated branch by the implicit 
function theorem. To prove the second part of the theorem, remark that N{H) keeps Vh invariant (as a 
set). Either the group N{H)/H = 1 or N{H)/H ~ Z2. This is because Q \s a, group of isometries whose 
action in V is orthogonal (orthogonal matrices). In the second case /(•, A) is an odd function which implies 
that the bifurcation is a pitchfork. D 

Of course it may also happen that solutions generically bifurcate in the principal stratum of a fixed point 
subspace which has dimension greater than 1. This point will be addressed in the example of the next 
subsection. 

Definition 3.3. Any generic T -periodic steady-state satisfying the hypotheses of the Equivaraint Branching 
Lemma is called an H-planform. 

Finally, if x' = g ■ x, then the isotropy subgroup oi x' is H' = g ■ H - g^^ . The image of x under Q is called its 
5-orbit. Equilibria or more generally flow-invariant sets belonging to the same group orbit share the same 
properties. 

3.2 Example: the octagonal lattice 
3.2.1 The octagonal lattice and its symmetries 

The octagonal lattice group F is generated by the following four hyperbolic translations (also called boosts), 
see [7]: 

\ ^2 + 2^/2 1 + ^/2/ 

and gj = ?'j7r/45o^-jir/4; i ~ Ij 2, 3, where r^ indicates the rotation of angle Lp around the origin in D. The 
fundamental domain of the lattice is a regular octagon C The opposite sides of the octagon are identified 
by periodicity, so that the corresponding quotient surface D/F is isomorphic to a " double doughnut" (genus 
two surface) [7]. 

We now determine what is the full symmetry group of the octagonal lattice, or equivalently, of the surface 
D/F. Clearly the symmetry group of the octagon itself is part of it. This is the dihedral group Dg generated 
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Figure 5: Tesselation of the regular hyperbolic octagon with the triangle T of Figure 2, colored in purple 
in the plot. We define two points S and S. S is the center of the rotation ct by tt (mod F), see text in 
subsection 3.2.1. S is the center of the rotation a hy tt (mod F), see text in subsection 3.2.1. 

by the rotation r.„-f4 and by the reflection k through the real axis, but there is more. We have seen in Section 
2.2 that the group G is generated by the reflecitons through the edges of the elmentary trinagle tiling D/F. 
The smallest triangle (up to symmetry) with these properties is the one shown in Figure 2. It has angles 
7r/8, 7r/2 and 7r/3 at vertices P ~ O (the center of D), Q, R respectively, and its area is, by Gauss-Bonnet 
formula, equal to 7r/24. There are exactly 96 copies of r filling the octagon, hence \Q\ — 96. Figure 5 shows 
this tessellation of O by triangles. 

The index 2 subgroup of orientation preserving transformations has 48 elements. In [9] it has been found 
that G ~ GL{2, 3), the group of invertible 2x2 matrices over the field Z3. In summary: 

Proposition 3.2. The full symmetry group Q o/D/F is G U nG where G ~ GL{2, 3) has 48 elements. 

The isomorphism between G'L(2,3) and G can be built as follows. We use the notation Z3 = {0, 1, 2} and 
we define: 

• p the rotation by 7r/4 centered at P (mod F), 

• a the rotation by tt centered at Q (mod F) 

• e the rotation by 27r/3 centered at R (mod F). 

Then we can proceed with the following identification. 



P = 



since these matrices satisfy the conditions p^ — a^ ~ e^ = Id and pae — Id. Note that p'^ = —Id where Id 
is the identity matrix. We shall subsequently use this notation. The group GL{2, 3) is described e.g. in [49]. 
The full symmetry group Q is generated by G and k, the reflection through the real axis in D. We further 
define 
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• k' = pK the reflection through the side PR of the triangle r, 

• k" = en' — an the reflection through the third side QR. 

The group Q and its representations have been studied with the help of the computer algebra software GAP 
[34]. Details are found in [18]. The main result is summarized in the next proposition and table. 

Proposition 3.3. There are 13 conjugacy classes in Q, hence 13 irreducible representations of Q . 4 of them 
have dimension 1, 2 have dimension 2, 4 have dimension 3 and 3 have dimension 4- Their characters are 
denoted Xj ■ j — 1, .-., 13. The character table is shown in table 1. Moreover all these representations are real 
absolutely irreducible. 



Class # 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


Representative 


Id 


P 


2 

p 


-Id 


a 


e 


— e 


K 


k' 


an 


pan 


en 


—en 


Xi 


1 


1 


1 


1 


1 






1 


1 


1 


1 






X2 


1 


-1 


1 


1 


-1 






1 


-1 


-1 


1 






X3 


1 


-1 


1 


1 


-1 






-1 


1 


1 


-1 


-1 


-1 


Xi 


1 


1 


1 


1 


1 






-1 


-1 


-1 


-1 


-1 


-1 


X5 


2 





2 


2 





-1 


-1 


-2 








-2 






X6 


2 





2 


2 





-1 


-1 


2 








2 


-1 


-1 


X7 


3 


1 


-1 


3 


-1 








-1 


-1 


1 


3 








X8 


3 


1 


-1 


3 


-1 








1 


1 


-1 


-3 








X9 


3 


-1 


-1 


3 


1 








1 


-1 


1 


-3 








Xio 


3 


-1 


-1 


3 


1 








-1 


1 


-1 


3 








Xii 


4 








-4 





-2 


2 




















Xl2 


4 








-4 





1 


-1 














%/3 


-73 


Xl3 


4 








-4 





1 


-1 














-V3 


Vi 



Table 1: Irreducible characters of Q. Notation: d = tat ^ (rotation by tt centered at S in Fig. 5) 

Proof. We only show here the second part of the proposition (see [18] for the first part). 
Absolute irreducibility is obvious for the one dimensional representations and it follows as a corollary form 
the next proposition 3.4 for representations X5 to xw- 

It remains to prove the result for the four dimensional representations xiii X12 and Xi3- For this we consider 
the action of the group Dg generated by p and k (the standard symmetry group of an octagon), as defined 
by either one of these 4D irreducible representations of Q. We observe from the character table that in 
all cases, the character of this action is x{p) — 0, x(p^) = 0, xi^Id) = —4, x{p^) = (p and p^ are 
conjugate in Q), and x('*) = x('*') — 0- We can determine the isotypic decomposition for this action of D^ 
from these character values. The character tables of the four one dimensional and three two dimensional 
irreducible representations of Dg, can be computed easily either by hand (see [57] for the method) or using a 
computer group algebra software like GAP. For all one dimensional characters the value at —Id is 1, while 
for all two dimensional characters, the value at —Id is —2. Since x(— /d) = —4, it is therefore not possible 
to have one dimensional representations in this isotypic decomposition. It must therefore be the sum of 
two representations of dimension 2. Moreover, since xip) — x(p^) = x(p'^) — Oj it can't be twice the same 
representation. In fact it must be the sum of the representations whose character values at p are ^2 and —\f2 
respectively. Now, these representations are absolutely irreducible (well-know fact which is straightforward 
to check), hence any Dg-equivariant matrix which commutes with this action decomposes into a direct sum 
of two scalar 2x2 matrices A/2 and pli where A and p are real. But the representation of Q is irreducible, 
hence A = /i, which proves that it is also absolutely irreducible. D 
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The 2 and 3 dimensional irreducible representations of Q must act as some particular irreducible represen- 
tations of finite subgroups of 0(2) (planar case) and of 0(3) (dimension 3). The next result specifies these 
actions. The idea is that these representations are not faithful. In other words for every j — 5, . . . , 10 
the kernel kerxj = f^j ¥" {0}- Let Xj be the representation induced by the projection G — >• Q/Hj. Then 
Xj is isomorphic to some irreducible representation of either a dihedral group D„ (symmetry group of the 
regular n-gon) in the 2-d case or the symmetry group of a platonic body in the 3-d case. We write O for the 
octahedral group (direct symmetries, i.e. rotations of a cube), and T^ for the full symmetry group of the 
tetrahedron. These two groups are isomorphic and have 24 elements. We also write Z2 for the 2-element 
group generated by the antipodal reflection in M"^. 

Proposition 3.4. (i) X5 is isomorphic to Dg acting in M^; 

(ii) xe is isomorphic to D3 acting in M? ; 

(Hi) X7 is isomorphic to O acting in M'^; 

(iv) xs is isomorphic <o O x Z2 acting in M.^ (full symmetry group of the cube); 

(v) X9 is isomorphic to the action ofTd x ^2 in ^^ ; 

(vi) Xw is isomorphic to the action ofTd in M"^. 

Proof. The proof of this proposition is given in [18]. We show a slightly different proof for the cases (iii) 
to (vi). We see from the character table that in X7 to XiO: the element —Id acts trivially. Let us write 
0*2 = {Id, —Id}. Then GL{2, 3)/02 is a 24 element group which is known to be isomorphic to the permutation 
group 5*4, hence to O. Therefore G/C2 ~ O x Z2. It follows that the four irreducible representations X7 to 
Xw must be isomorphic to the four irreducible representations of dimension three of the symmetry group of 
the cube O x Z2 . Comparison of the character tables of these representations leads to the result (one can 
make use of GAP to obtain these tables) . D 

The actions listed in the proposition are absolutely irreducible. It follows that the representations X5 to Xio 
are also absolutely irreducible. 

3.2.2 Classification of H-planforms and bifurcation diagrams 

It follows from Proposition 3.4 that the generic bifurcation diagrams for all representations Xi to Xio are 

classical. These diagrams are well-known for D3 and Dg, [35] and they are also known for the 3-dimcnsional 

cases. 

Let us for example consider cases (iii) and (iv). By Theorem 3.1 each axis of symmetry of O and O x Z2 

gives rise to a branch of equilibria with that symmetry. These are the symmetry axes of a cube and they here 

of three types: 4 fold symmetry (passing through the centers of opposite faces), 3 fold symmetry (passing 

through opposite vertices) and 2 fold symmetry (through opposite edges). Moreover each symmetry axis is 

non trivially mapped to itself by some transformation in the group. It follows that there are three different 

types of pitchfork branches of equilibria in these cases. For the stability and (generic) non existence of other 

branches we refer to [56]. 

What remains to do in these cases is to actually compute the bifurcated states, which is quite involved since 

the eigenf unctions of Ap/r can't be expressed in an explicit manner. This part will be discussed in the next 

subsection. 

In the rest of this paragraph we study the case of 4 dimensional representations xii to X13, which do not 
belong to the list in Proposition 3.4. This was studied in details in [30], to which we refer for proofs. Here 
we outline the main results and sketch the proofs. 
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First note that the dimension of the fixed point subspaces of isotropy subgroups for a representation Xj can 
be determined from the character table 1 thanks to the formula [15, 35] 



dim(0 = ^^X,(/^) 



(3.5) 



h£H 



In order to find the H-planforms we therefore need to find those isotropy subgroups such that -mr X^/ig// Xj (^) = 
1. However it is important also to determine all the other (classes of) isotropy subgroups because (i) branches 
in higher dimensional fixed point subspaces should not be a priori excluded and (ii) this is useful for the 
study of the stability conditions and local dynamics. In [18] all subgroups of Q have been determined (using 
GAP), a necessary step to compute the istropy subgroups using formula 3.5. The list is cumbersome, so we 
simply list the isotropy subgroups for the representations xn to Xi3 and refer to [18] for the details. We 
shall use the notation a = p^ap'^^, which is the rotation by tt around point S in Figure 5. 



Definition 3.4. 



63 - 



(cr, k) = {Id, (7, K, k"} 
(a, k) = {Id, a, K, —p^n"p^^{ 
(e, k!) — {Id, e, e^, k', ek's^ , e^n'e} 
{<T,e) = {/d, e, e^,cr, eae^, e^CTe} 



The meaning of these groups for the octagonal lattice is not straightforward and will be explained in the 
next subsection. 

1. The xn o.'^'i- Xi3 cases. 

The (conjugacy classes of) isotropy subgroups are ordered by set inclusion and we call lattice of isotropy 
types the corresponding graph. 

Lemma 3.1. The lattices of isotropy types for the representations X12 o.'^^d Xiz '^'"e identical and are shown 
in Figure 6. The numbers in parentheses indicate the dimension of corresponding fixed-point subspaces. 

(0) 



(1) 



(2) < K > 



(4) 




Figure 6: The lattice of isotropy types for the representations X12 and X13. 



Remark 3.1. (i) The difference between X12 and X13 is quite subtle as we see from the character table and 
their geometrical propreties are similar as well as their bifurcation diagrams for Equation (3.2). 
(ii) We can see on the lattice of isotropy types that each fixed-point plane contains at least 2 or 3 axes 
of symmetry. The exact number is 4- For example Fix{{a)) contains one copy of Fix{C2K), one copy of 
Fix{C2f^) and two copies of Fix{Ti^). 

Thanks to Theorem 3.1 each axis of symmetry contains bifurcated solutions, and moreover the branches 
are pitchfork because one can easily check that for each such isotropy subgroup H (as defined in 3.4), 
N{H)/II ~ Z2. In order to study the stability of these solutions and to look for other equilibria we need 
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to know the Taylor expansion of the bifurcation equation (3.2) up to some order large enough to remove 
degeneracies. It turns out that in this case cubic order in X is enough. However here again the calculations 
are cumbersome and we refer the reader to [30] for details. 

We first need to choose suitable coordinates for X e M* and we do so for the representation Xi2 (we know 
similar results hold for xi3- Let {zi,zi, Z2, Z2) be the coordinates which diagonalize the matrix of the 8-fold 
symmetry p in M''. The equation (3.2) truncated at order 3 in these coordinates reads 

ii = [A + a(|zi|2 + |z2n]zi+fe[V3(3z2 + £|)zi-i(z|+3z2)z2] (3.6) 

i2 = [A + a(|zi|2 + IzsH] Z2 + b\y3 (3z| + z^) z^ + i (z? + 3Z2') 21] (3.7) 

and their complex conjugates, where a, h are real coefficients. These equations can be written as a gradient 

system. 

We can now state the main result for representations xvi ^^id Xi3- 

Theorem 3.2. Provided that (a, b) £ V ^ {{a,b) e R^ \ 3a + 2bVS < and 3a + lObVS < 0}, the following 
holds for Equations (3.6)-(3.7). 

(i) No solution bifurcates in the principal strata of the planes of symmetry. 

(ii) The branches of equilibria with maximal isotropy are pitchfork and supercritical. 

(Hi) If b > (resp. b < 0), the equilibria with isotropy type C^^i (resp. D3J are stable in M'*. Branches 
with isotropy C2K; and C2„ are always saddles. 

We didn't prove the non existence of bifurcated equilibria with trivial isotropy, however there is a strong 
evidence that such solutions are generically forbidden for equations (3.6) and (3.7). 
We now turn to the xn case. 

2. The xii case. 

Lemma 3.2. The lattice of isotropy types for the representation Xu is shown in Figure 7. The numbers in 
parentheses indicate the dimension of corresponding fixed-point subspaces. 

G 



(1) r r' 



(2) < K > < a > < K > 



Id 

Figure 7: The lattice of isotropy types for the representation xii- 

We see that this representation is somewhat different from the other 4 dimensional representations of Q, and 
indeed it leads to quite different bifurcation diagrams. In particular il allows for the bifurcation of equilibria 
with isotropy a or k! under generic conditions on the coefficients of terms of order 5 in the Taylor expansion 
of the bifurcation equation. Indeed there is no term with even order and the only term of order 3 is the 
"radial" one: ||X|pX, which can't provide pattern selection. The next order is 5 and one can show that 
there are 4 independant ^-equivariant terms at this order. Details are provided in [30]. Among these terms, 
three are gradients while one is non gradient. It was shown in [30] that when the latter is non zero, a robust 
heteroclinic cycle can bifurcate from the trivial state. However this phenomenon, which is interesting in the 
context of non gradient systems, does not occur for the Swift-Hohcnberg equation. 
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As before let (zi, zi, Z2, Z2) be the coordinates which diagonahze the matrix of the 8-fold symmetry p in M'*. 
The equation (3.2) truncated at order 5 and with no non gradient terms reads in these coordinates 



d 

Zi = [X + A{\zi\^ + \Z2\^)] Zi + —Q{zi,Zi,Z2,Z2) 

d 

h = [X + A{\zi\^ + \Z2\^)] Z2 + -^Q{zi,Zi,Z2,Z2) 



(3.8) 
(3.9) 



where Q is defined as follows (C.C. means "complex conjugate"): 

Q{zi,zi,Z2, Z2) = a{\zi p + \z2ff + b{zfziZ2 + 2zizfz^ - ziz^Z2 - 2zfz2z1 + C.C.) + d{ziz^ - zjza + C.C). 

Compared with the equations in [30], we see that a term with coefhcient c in the latter does not appear 
here. This is due to the gradient structure and the correspondance between the two equations is obtained 
by setting c = — 26 in [30]. 

In these coordinates the equations of the fixed-point subspaces are 

- Fix(CT): Z2^iil + V2)zi, 

- Fix(K): zi — izi and Z2 — ~*^2; 

- Fix(K'): zi = ^(?; - l)zi and 23 = ^(i + 1)^2, 

and the intersections of these planes are the fixed point axes for the isotropy subgroups C2K and C'21^. Using 
these informations one can show the following 

Theorem 3.3. Provided that A < in (3.8)-(3.9), the branches of equilibria with isotropy C2K and C2K are 
pitchfork and supercritical. Moreover: 

(i) Fix^K) contains two copies of each type of axes of symmetry and there is generically no other solutions 
bifurcating in this plane. 

(ii) Fix{(j) and Fix(K') contain each only one copy of both types axes of symmetry. Moreover in either 
plane solutions can bifurcate in the principal stratum, depending on the values of the coefficients b and 
d in the bifurcation equations (3.8)-(3.9). The bifurcation and stability diagrams in terms of b and d 
are shown in Figure 9. The regions TH-III which are referred to in this figure are defined in Figure 8. 
Region IV is similar to II (bifurcation of solutions in the principal stratum). 
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(a) In Fix(cr), solutions with isotropy cr exist in 
the regions // and IV. We have set x = —b 
and y = d. 



(b) In Fix(K'), solutions with isotropy k' exist 
in the regions // and IV . We have set x = ~2b 
and y = h + d. 



Figure 8: Conditions of existence of solutions with isotropy a and n' (in blue). 



19 






C2k 



HI 




^ C2k 



Figure 9: Changes of phase diagram in Fix(cr) and Fix(K') as the coefhcients {b, d) pass from regions / — 
// — /// in Figure 8(a) (for Fix(cr)) and in Figure 8(b) (for Fix(K')). Sj^ and S2 indicate solutions with 
isotropy a or k' . The case II / III corresponds to coefficient values at the boundary between region // and 
/// ( saddle- node bifurcation of equilibria with isotropy a or k'). 
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3.2.3 Computation of H-planforms 

It follows from the definition that H-planforms are eigenf unctions of the Laplace-Beltrami operator in D 
which satisfy certain isotropy conditions: (i) being invariant under a lattice group F and (ii) being invariant 
under the action of an isotropy subgroup of the symmetry group of the fundamental domain D/F (mod F) . 
Therefore in order to exhibit H-planforms, we need first to compute eigenvalues and eigenfunctions of Ad in 
D, and second to find those eigenfunctions which satisfy the desired isotropy conditions. In this subsection, 
we use the notations /i and ^ for eigenvalue and eigenfunction of Ad: 

-Ad*(z) = A**(z), Vz e O. 

Over the past decades, computing the eigenmodes of the Laplace-Beltrami operator on compact manifolds has 
received much interest from physicists. The main applications are certainly in quantum chaos [4, 5, 7, 23, 61] 
and in cosmology [22, 40, 50]. 

In order to find these H-planforms, we use the finite-element method with periodic boundary conditions. 
This choice is dictated by the fact that this method will allow us to compute all the first n eigenmodes and 
among all these we will identify those which correspond to a given isotropy group. 
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Figure 10: Plot of the eigenfunction of the Laplace-Beltrami operator in the Poincare disk D with tj-symmetry 
with eigenvalue fx = 23.0790. 

As there exists an extensive literature on the finite element methods (see for an overview [1, 19]) and as 
numerical analysis is not the main goal of this review, we do not detail the method itself but rather focus on 
the way to actually compute the eigenmodes of the Laplace-Beltrami operator. We mesh the full octagon 
with 3641 nodes in such a way that the resulting mesh enjoys a Dg-symmetry We implement, in the finite 
element method of order 1, the periodic boundary conditions of the eigenproblem and obtain the first 100 
eigenvalues of the octagon. Our results, as reported in [18], are in agreement with those of Aurich and 
Steiner [4]. 

We plot in Figure 10, the first eigenfunctions of the Laplcae-Beltrami operator with full octagonal symmetry Q 
with non zero eigenvalue. It is associated to representation xi- Next, in Figure 11, we plot the corresponding 
eigenfunctions of the Laplace-Beltrami operator associated to the lowest non-negative eigenvalue ^ = 3.8432 
with multiplicity 3. We identify each solution by its symmetry group. It is clear that 11(b) and 11(c) can be 
obtained from 11(a) by hyperbolic transformations. From the definitions of go in (3.4), we see that g^ — Uro 
with ro =\n(l + V2 + ^2 + \/2 ) . If we define 7fc G g by: 

7fc = rotfc^/4ar„/2rot_fc^/4 



(3.10) 
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(a) H-planform \I/3 with symmetry Dg 





(b) H-planform 'I'l with symmetry 71 ■ Dg 



(c) H-planform \I'2 with symmetry 72 ■ Dg. 



Figure 11: Plot of the eigenfunctions of the Laplace-Beltrami operator m the Pomcare disk D associated to 
the lowest non-negative eigenvalue ^ = 3.8432 corresponding to the irreducible representation xs- In (a) we 
also plot the octagon (black line) and in (b),(c) its image by 71,72 respectively. 
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then Figure 11(b) (resp. 11(c)) is obtained from 11(a) by applying 71 (resp. 72). A gallery of eigenfunctions 
can be found in [18, 30]. 

3.2.4 Case study: computation of the bifurcation equations for irreducible representation xs 

For most of the irreducible representations, it is practically impossible to compute the reduced equation 
(3.2) given by the center manifold theorem. The main reason being that we only know numerically the 
eigenfunctions of the Laplace-Beltrami operator. It turns out that we have been able to successfully conduct 
this computation in the case of 3 dimensional irreducible representation xsj when the evolution equation is a 
neural field equation [29] . In that particular study, the choice of xs was dictated by the simple interpretation 
of H-planforms given in Figure 11 in terms of preferred orientations within a hypercolumn of the visual cortex 
and allows the computation of geometric visual hallucinations across the cortex [29] . As already explained in 
the previous section, the lowest non-negative eigenvalue /i — 3.8432 is of multiplicity 3 and associated to the 
irreducible representation xg. We then select the parameter a in equation (1.1) to be equal to ac such that 
for iXf. — n — 3.8432 we have a^ = ^c- Restricting to the class of F-periodic functions, this ensures that the 
first modes which bifurcate from A = are associated to the irreducible representation xsj all other modes 
being damped to zero (see discussion in section 3.1.1). We then rewrite the Swift-Hohenberg equation (1.1) 
with a = ac 

ut = -{al + l^ofu + \u + N{u). (3.11) 

Proposition 3.5. For the three dimensional irreducible representation xs of Q , the isotropy subgroups with 
one dimensional fixed point subspace are the following: 

Ds = {p, k) 

D2k = {-Id,a,K). 

In Figure 12, we represented the different axes of symmetry of the cube with isotropy subgroups given in 
proposition 3.5. Note that planforms in Figure 11 correspond to the three coordinate axes of the cube. 



(0,C3, -C3) 



71 -D, 




(0,0, C3) 



72 -Ds 



(C3,-C3,C3) 



a'^^ 
/- 



C2 



Cl 



Figure 12: Maximal isotropy subgroups Dg, CgK' and D2K of O k Z2. The axes 71 • Dg and 72 • Dg are copies 
of Dg by the elements 71,72 € G (see (3.10)). The plane (0,02,03) (resp. (01,0,03)) has symmetry D2 (resp. 
D'). 
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If we denote ^i the H-planform in Figure 11(b), ^2 the H-planform in Figure 11(c) and ^3 the H-planform 
corresponding to the symmetry group Dg in Figure 11(a), then (^1,^2,^3) is a basis for the irreducible 
representation xs- This can be easily seen through the identification of each H-planform to the three 
coordinate axes of the cube in Figure 12. Note that we have normalized planforms such that: 

(^^>^i) = ^^^^(^)*,(2)dm(z) = S.,,. 

We rewrite equation (3.11) as a dynamical system in the infinite-dimensional phase space X — L (O) 
consisting F-periodic functions: 

Ut = Lam +N[u) 

with 

The linear part L^ is a closed linear operator with dense domain y = 'i{t^^{0) where 

B^%AO) = {ue HL(D) I «(7 • z) = u(z), V(7, z) G g X D} 

and the nonlinear map A/" : 3^ — > 3^ is C for all positive integer p. The linear operator Lq is closed in X , 
with dense and compactly embedded domain y. Then Lq has a compact resolvent, and its spectrum cr(Lo) 
is purely point spectrum, only. Moreover we have cr(Lo) n iM = {0}. The eigenvalue is geometrically and 
algebraically triple (consequence of the absolute irreducibility of the representation of Q in that eigenspace) . 
We denote V the corresponding three-dimensional eigenspace. A direct calculation also shows that there 
exists Wo > and C > such that for all |a;| > wq, we have the resolvent estimate: 

\\{iu-U)-^\\c(x)<^,- 

We can apply the equivariant center manifold reduction introduced in 3.1.2. So there exists a map $ : 
V^ X M — > V^ such that all solutions u that lie on the center manifold can be written 

u = xi'^i -t-a;2*2 + 2:3*3 -I- $(a;i,a;2,a;3. A), 

and (ci, C2, C3) satisfies the following reduced system: 

\xi + [a{x\ + x\) + bxf\ xi + h.o.t. 



dt 

-^ = \x2+[a{xl + xl) + bxl]x2+\^.o.t. (3.12) 

— — = Ax3 + [a{x\ + xl) + bxf\ 2:3 + h.o.t. 



We refer to [56] for the computation of the normal form (3.12) and for a review on bifurcation problems with 
octahedral symmetry. 

Taylor expanding the map $: 

$(a;i,a;2,X3,A) = ^ x^a:^x[jA™$^si^ -H h.o.t 

0<r+s+;+r?j.<3 

and denoting TZ: 

TZ{u, A) = Am +J\f{u) — Tlii{u, A) + TZ2q(u, u) + 7?.3o(u,m,m) 

with 

TZii{u, A) = Xu 
TZ2o{u,v) — vuv 
TZ30o{u,v,w) — —rjuvw, 
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37^3o(*l,*3,*3),*l) 
.*3). 



(3.13) 



we obtain the following system of equations: 

= -Lo$0020 - '7^2o(*3, ^3) 

= -Lo$ioio - 27e2o(*i, ^3) 

a = (27^2o(«'oo20, *i) + 27^2o($lolo, ^3) 

b = (27^2o(*3, *0020) + 7^30(*3, *3, ^3) 

Here (•, •) is the scalar product on Lpg^(C'). 

In order to solve the two first equations of the previous system, we need to know if the functions '^3{z)'^3{z) 
and ^1/1(2)^1/3(2;) can be expressed as a linear combination of eigenfunctions of the Laplace-Beltrami operator 
on C In general, it is very difficult to obtain these expressions because the eigenfunctions are only known 
numerically and one needs the computation of the associated Clebsch-Gordan coefficients. It turns out that 
in our case we have been able to conjecture and numerically verify the following relations: 



^i(z)*3(z) 



1 



V3 



^^k{z) 



*lW = ^*?e^"W + l 



where the corresponding isotropy subgroups are given by: 



Dk=< 



-Id, p K, p a > and Dg^ =< P, P crp 



,K > . 



.Di 



The notation VP^^j^o" (2) means that the product '^i{z)'i'3{z) is an eigenfunction of the Laplace-Beltrami op- 
erator associated to the irreducible representation xio with isotropy subgroup D2K. Similarly the notation 
(z) stands for an eigenfunction of the Laplace-Beltrami operator associated to the irreducible represen- 



X6 



D' 



tation xe with isotropy subgroup Dg^. Furthermore we have normalized '^xw ^^'^ ^^r**" such that: 



X6 

\^Xio : *Xio / \ xe ' Xe ' 

In Figure 13, we plot the eigenfunctions '^xw '"^^^ ^P**" of the Laplace-Beltrami operator in the octagon 
O. One interesting remark is that the product ^1^3 corresponding to the three dimensional irreducible 
representation xs produces an eigenfunction associated to another three dimensional irreducible representa- 
tion: xw whereas ^| is the linear combination of the constant function which has G as isotropy subgroup 
and thus corresponds to xi and the eigenfunction ^^g^" which is associated to two dimensional irreducible 

D' 



representation xe- We denote Hxw a-^^d mP*" the corresponding eigenvalues 



-Ad*°4» = fixk'Sxk and - AdvI'^- = m°-* 



Xe 



With these notations, the two first equations of system (3.13) give 
$0020 = Span(*i,*2,*3) + i^ 
$1010 = Span(*i,*2,*3) 



* 



2iy 



A straightforward but lengthly calculation gives the expression of the coefficients a and b in the reduced 
equation (3.12) 

■" 2 8 4 



L"c 5ial-p^ty Hal~p'^,l-y 



V 



b= '' " ^^ 

< 2b{al-p^tY\ 

From the analysis derived in [56] , we have the following result 



61 
25' 



(3.14) 
(3.15) 
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(a) Plot of *°4" • 



(b) Plot of ^^g*" . 



Figure 13: Plot of the eigenfunctions of the Laplace-Behrami operator in the octagon O corresponding to the 
irreducible representations xio with eigenvalue /ixi^o" = 15.0518 (left) and xe with eigenvalue A^xe**" ~ 8.2501 
(right). 

Theorem 3.4. The stability of the branches of solutions corresponding to the three maximal isotropy sub- 
groups given in proposition 3.5 is: 

(i) the Dg branch is stable if and only if a < b < 0, 
(a) the CgK' branch is stable if and only if 2a + b < and b ~ a < 0, 
(Hi) the D2k branch is never stable. 

For each value of {i',r]) in [0,5] x [0,5] we have numerically computed the coefficients a and b given in 
equations (3.14)-(3.15) and then checked if the stability conditions in Lemma 3.4 are satisfied. Our results 
are ploted in Figure 14. We can see the different regions of the plane (z^, r]) where the branches Dg and Cere' 
are stable: in bright gray the region where Dg is stable and in dark gray the region where CgK' is stable. 
Note that there is a whole region of parameter space where all the branches are unstable. 



4 Radially localized solutions 

Recently, there has been much progress made in understanding radially localized solutions in the planar 
Swift-Hohenberg equation. For the Swift-Hohenberg equation near the Turing instability, three types of 
small amplitude radially symmetric localized solutions have been proven to exist: a localized ring decaying 
to almost zero at the core, a spot with a maximum at the origin (called spot A) and a spot with minimum at 
the origin (called spot B); see [51, 53-55]. The proofs rely on matching, at order 0(l/r^), the "core" manifold 
that describes solutions that remain bounded near r = with the "far-field" manifold that describes how 
solutions decay to the trivial state for large r. The core manifold is found by carrying out an asymptotic 
expansion involving Bessel functions while the far-field manifold is calculated by carrying out a radial normal 
form expansion near ?' = oo. The main technical difficulty is that the far- field normal form is only valid 
down to order 0{l/r) and so the manifold has to be carefully followed up to order 0(l/r^). Localized rings 
occur due to a localized pulse in the far-field normal form equations and require that the bifurcation of rolls 
at r = oo is subcritical. Spot A solutions occur due to the unfolding of a quadratic tangency of the core 
manifold and the cubic tangency of the far- field manifold with the trivial state at onset; see [51, Figure 4]. 
The spot B state is formed by 'gluing' the spot A and localized ring solution. Crucially, all these localized 
radial states are L -functions that can not be found via a Lyapunov- Schmidt or center manifold reduction. 
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V 




Figure 14: Regions of the plane {v, rf) where the branches Dg and Cbk' are stable: in bright gray the region 
where Dg is stable and in dark gray the region where CgK' is stable. 

In this section, we will only be interested in the existence of spot A type of solutions for the Swift-Hohcnberg 
equation (1.1). The proof of the existence of such solutions will closely follow the one presented by Faye et 
al. [33] for neural field equations on the Poincare disk. In the hyperbolic case, the major difficulty comes 
from the fact that it is not clear how to define the core and far-field manifolds in order to carry out the 
matching. It turns out that the far-field manifold is easier to define than in the Euclidean case since there is 
no bifurcation in the far-field. However, calculating the core manifold is significantly more involved than in 
the Euclidean case and constitutes the main challenge in the existence proof of spots in hyperbolic geometry. 
From now on we shall use the terms, spot and bump interchangeably to refer to the spot A states. 

4.1 Notations and definitions 

Throughtout this section, we work in geodesic polar coordinates z = {T,ip) G D, with z — tanh(r/2)e*'^. In 
these coordinates, the measure element defined in equation (2.2) is transformed into dm(z) — smh{T)dTdip. 
Furthermore, in order to fix ideas, we set the value of a in equation (1.1) to be equal to 1: 

ut^ -{1 + ADfu + Xu + uu^ -r]u^. (4.1) 

The Laplace-Beltrami operator defined in equation (2.5) can be written in geodesic polar coordinates as 

(4.2) 



Ad = ^ + coth(T)— + sinh(T)-2_. 



We define the radial part of the Laplace-Beltrami operator to be 

A..^+coth(r)-. 



(4.3) 



Stationary radial solutions u{t) of equation (4.1) depend only on the radial variable r = do{z, 0) and therefore 
satisfy the ordinary differential equation 







(1 + Ar)^u + Xu + mi^ - T^U^ 



(4.4) 



We are interested in finding localized solutions u(r) of (4.4) that decay to zero as t — ?► cx) and that belong 
to the functional space L^(IR+,sinh(T)dr). We shall therefore seek for such solutions for A < 0, where the 
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background state u = is stable with respect to perturbations of the form e°'*ep^b{z) (see 3.1.1). In fact, our 
results are restricted to < — A ^ 1, and we shall construct localized radial solutions with small amplitude 
that bifurcate from u = at A = into the region A < 0. 

In the Euclidean case, Bessel functions play a key role in this analysis of equation (4.4) close to the core 
T = 0. In the hyperbolic setting, the analog of the Bessel functions are the associated Legendre functions of 
the first king V^ and second kind Q^. We now recall their definition (see Erdelyi [27]). 

Definition 4.1. We denote V^lz) and Q^{z) the two linearly independent solutions of the equation 

£■ d / B'^ \ 

(1 - z^)-—u{z) ~ 2z—u(z) + a(a + 1) - -^—^ u{z) = 0. (4.5) 

dz^ dz \ 1 — z^ / 

"Pf (z) and Qa{z) are respectively called associated Legendre function of the first and second kind. For f3 — 0, 
we use the simplified notation Va{z) = V'^{z) and Qa{z) = Qq(z). 

4.2 Main result 

We can now state the main result of this section. 

Theorem 4.1 (Existence of spot solutions). Fix v ^ Q and any ry G M, then there exists A.^ < such that 
the Swift-Hohenberg equation (4.1) has a stationary localized radial solution u{t) G iy^(M+, sinh(r)(iT) for 
each A G (A», 0).' these solutions stay close to u = and, for each fixed t» > 0, we have the asymptotics 



w(t) = cV|A| 7'_i^.y3(coshT) + 0(A) as A^O, (4.6) 

2 ' ^ 2 

uniformly m < r < r* for an appropriate constant c with sign(c) — sign(z/). 

As for the Euclidean case, this theorem states that spots bifurcate for any value of i^ 7^ 0, regardless of the 
value of 77. The matching arguments in 4.3.2 yield a similar theorem for the case A > where one finds 
a bifurcating branch of solutions also given by (4.6). However, these solutions are not _L^(]R+, sinh(r)dr)- 
f unctions. 

Note that spots are initially unstable with respect to the PDE dynamics of the Swift-Hohenberg equation 
(1.1). However, our numerical continuation results indicate that spots stabilize in a certain region in the 
(A, i^)-parameter space, at least when the nonlinear part of equation (1.1) is Af{u) = vu^ — u^ . 

In the following subsection we sketch the proof of Theorem 4.1. It closely follows the one given in [33] to 
which we refer for further details. 

4.3 Proof of Theorem 4.1 
4.3.1 The equation near the core 

We first rewrite equation (4.4) as a four dimensional system of non-autonomous differential equations to 
yield 

drUi = M3, 
9rU2 = U4, 

(pi -t- coth(r)9T- -ir\)ui— U2, 

(pi -t- coth(r)9T 4- 1) U2 = Aui + vu\ — rju\, 

and we may rewrite (4.4) as a spatial dynamical system 

Ur=A(T)U + F{U,\), (4.7) 
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A{r) = 



, :f{u,\) = 



\ 



\\ui + vu\ — riu\/ 



ith 

/ 1 ^ \ 

1 

-1 1 -coth(r) 

V -1 -coth(r)/ 

and U — {ui,U2,us,U4)'^ . 

We begin our analysis by characterizing all small radial solutions of (4.7) that are bounded and smooth 
in the interval [0,ro] for any fixed tq. We set A = and linearize (4.7) about U — to get the linear 
system drll = A{t)U. Proposition 4.f defines the four linearly independent solutions of the linear system 
drU = A{t)U. 

Proposition 4.1. The linear system drU = A{t)U has four linearly independent solutions given by 

Viir) = (T'^^ (cosh t),0,7'1„ (cosh t),0)'^ 

V2{r) = {V2\T),Vao{coshT),Vi{r),Vi^{coshT)f 

Vsir) = (Q„„ (cosh T),0,Qi„ (cosh t),0)'^ 

Viir) = {Vl{T),Q^,{coshT),Viir),Ql^{coshr)f 

ao = -2+«^ (4.8) 



wher 



and 



V2{t) —Vaoicoshr) / 7'Q(j(coshs)QaQ(coshs)sinh(s)cfs — Q„Q(coshT) / (^'^^(coshs)) sinh(s)(is, 
Jo Jo 

V2 {t) ='Pa^{cos'\iT) I "Pao (cosh ,s)Qq,j, (cosh s) sinh(s)ds - Qi!j|^ (cosh t) / (^^^(coshs))^ sinh(s)ds, 
Jo Jo 

T r 2 r 

^ (■'') ~^Qo(coshr) / (^^^(coshs)) sinh(s)ds — Qcto(coshr) / T'qq (cosh s)Qq(j (cosh s) sinh(s)(is, 
Jo Jo 

VI{t) ='Pa^^{cos'\iT) / (QQ„(coshs))^sinh(s)ds - Qi!,^(coshr) / T'og (cosh 3)2^^ (cosh s) sinh(s)ds. 
Jo Jo 

Proof. From equation (4.5), the associated Legendre functions ^'^(cosh-) and QQ,(cosh-) form a basis of 
solutions for the equation 

a^*(r) + coth(T)a^*(r) - a{a + l)*(r) = 0. 

If a = —\ +* 2 ^^cn we have a{a + 1) = —1. This implies that T^i ^5 (coshr) and Q_i v5(coshT) are 



solutions of 

From now ona = ao = — 5+* 
the equivalent system 



dl^{T) + coth(r)9^*(r) + *(t) = 0. 
^. The solution of linear system drU = A{t)U can be found by inspecting 

{d^ + coth{T)dr + 1) ui = U2, {dl + coth(r)ar + l) U2 = 0. 

Consequently, solutions Vi (r) — (Pq^ (cosh t),0,7'q^ (cosh r), O) andV3(T) = (Qq„ (cosh t), 0, Q^^ (cosh t), O) 
are found to be two linearly independent solutions, where we have used the relations 

dr (7'a,(coshr)) = 7'^„(coshr), dr (Qao(coshr)) = Qi^(coshr). 

In order to find the other two linearly independent solutions, we have to solve the two equations 

{dl + cofhrdr + l) u{t) = 'Pao(coshT), (4.9) 

(a^ + cothTa^ + 1) u{t) = Q„„(coshT). (4.10) 
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From formula in [27, page 123] 

pQ (cosh r)Q;!j, (cosh r) — T'^ (cosh r)Qc« (cosh r) = — . " , r > 0, v e 

and straightforward computations, we obtain that the solutions of (4.9) are given by 

u{t) = ciVao (cosh r) + C2 Qao (cosh r) + Pqq (cosh r) / T'ao (cosh s) Qag (cosh s) sinh(s)ds 

Jo 

r 2 

— QQo(coshr) / (7'Qg(coshs)) sinh(s)ds, 
Jo 

with ci,C2 € M two real constants. Equivalently, the solutions of (4.10) are given by 

m(t) = CsVaoicoshr) + C4Qao(coshT) +'Pao(coshT) / (QQ,Q(coshs))^sinh(s)fis 

Jo 

— QQo(coshr) / "Pqj, (cosh s)Qc«o (cosh s)sinh(s)ds, 
Jo 

with constants C3, C4 e M. 

Chosing only linearly independent solutions, we finally obtain the result stated in the proposition. D 

In Table 2, we summarize the expansions of the associated Legendre functions in the limits r ^> and 
T — > 00; Proposition 4.2 defines the constants in the r — > 00 limit. Thus, Vi{t) and V2(t) stay bounded as 
T —> 0, while the norms of V3(r) and V4(r) behave like Inr as r — ?► 0. We expect that the set of solutions 
of (4.7) that are bounded as r —> forms a two-dimensional manifold in M'* for each fixed r > 0. We 
denote the projection P™(to) onto the space spanned by Vi(to), V2(to) with null space given by the span of 
^3(to),V'4(to). 





r^O 


T — > 00 


^a„(cOsh.) 


1 + 0(t2) 


Co cos(4^ + $o)e-5 + 0(e-¥) 


^io(cOSh-) 


r(-i+0(r2)) 


Co cos(4^ + $0 + f )e-5 + 0(e-^) 


Qao (cosh-) 


(-l + 0(T2))lnr + 0(l) 


CiCos(^^ $1)6-5+0(6-^") 


Sio(cosh-) 


(l + 0(T2))rlnT-i+0(l) 


Ci cos(4^ - $1 + f )e-5 + 0(e-^) 



Table 2: Expansions of associated Legendre functions 7'^^(cosh-) and QQo(cosh-) for r —J- and r — >■ 00; see 
[27, 66]. Co, Ci, $0 and $1 are all real constants given in Proposition 4.2. 

Proposition 4.2. The constants Co,Ci,<I>o a?irf ^1, given in Table 2, are 



Co = 2, 



\ 



2\/3 



Stt tanh 



(4^) 



anrf Ci 



\ 



27r a/3 tanh ('^') 



$0 = arg 



T{^f' 



r k + ^' 



(i 



and <I>o + $1 = - ' 



Proof. See [33]. 



(4.11) 



(4.12) 



D 



We are now able to present the hyperbolic equivalent of Lemma 1 [51] for the center-unstable manifold 
^^(A), the set of bounded and continuous solutions of (4.7) close to r = 0. That is, fix tq > and Sq > 0: 

n;™(A) = {[/ e C°([0,ro],M'*) solution of (4.7) | sup \U{t)\ < Sq for |A| < ,5o}. 

0<T<To 
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Lemma 4.1. Fix tq > 0, then there exist constants 60, Si such that the set yV™(A) of solutions U{t) of 
(4.7) for which suPq<^<^|.j |t^('^)l < '^o is, for |A| < 5{), is a smooth two-dimensional manifold. Furthermore, 
U e yVf:"(A) with |Pf"(To)L/(To)| < Si if and only if 

U{to) = diViiro) + d2V2{To) + l/3(To)Oro(|A||J| + |dT) 

/ . ... X (4.13) 

+ V^{to) (- [1 + 0(1)] K + 0.„(|A||d| + |d2p + |di|3)) , 

with 

f°° 3 

X= / (7^Q,Q(coshs)) sinhs(is<oo, 

Jo 
for some d — {di,d2) G K^ with \d\ < Si, where the right hand side in (4.13) depends smoothly on (rf, A). 

Proof. We observe that four independent solutions to the adjoint problem drU ~ —A^{t)U are given by 
(see [33]): 

Wiir) = smhr{Qi^icoshT),Wi{T),~Qa.^,{coshT),Wt{T)f, 

W2iT) = sinh T(0,Qi„ (cosh r),0,-Qc«o (cosh r))"^, 

Wsir) = smhT{vl{coshT),W^{T),~Va.o{co8hT),W^{T)f 



Wiir) = sinh t(0,7'^^ (cosh r),0,-7'a„ (cosh r))"^^ 



where 



W?ir)^V^coshT) f V^coshs)Ql^{coshs)smh{s)ds - Ql^icoshr) [ {V^cosh s)f sinh{s)ds , 
Jo Jo 

Wi{t) =- Qaoicoshr) -Vaoicoshr) / T'^^j (cosh s)Q;!j^ (cosh s) sinh(s)ds + Qqq (cosh t) / (7'^^(coshs)) sinh(s)ds, 

Jo Jo 

Wiir) ^Vi^icoshr) [ {Qi^icoshs)f smh{s)ds - Qi^icoshr) [ 7'^„(coshs)Qi Jcoshs) sinh(s)ds, 
Jo Jo 

W^{t) ^-'Pao {cosh t) -Vao {cosh t) / (Q;!j^(coshs)) sinh(s)(is - ^^^(coshT) / T'^j^ (cosh s)^;!^^, (cosh s) sinh(s)ds. 

Jo Jo 

It foUows from 

sinhr (P^^ (cosh t)Qqo (cosh r) - Pq^ (cosh r)Q;^^ (cosh r)) = 1, 

that 

(1^,(t), Wj(r ))r4 =S^j i,j = 1,...,4, 

is independent of r. For a given d — (di, ^2) G M^, we consider the fixed-point equation: 

C/(r) = ^rf>,(r) + ^V^,(r) f {W,{s),HU{s),\))ds + Y.^,{t) f\w,{s),T{U{s),X))ds 

= ^djVjir) + ^Vj{t) Wj,4{s)MU{s),X)ds + ^Vj{T) WjA{s)MU{s),\)ds 

j=l j = l Jro j=3 "'0 

(4.14) 

onC"([0,To],M''), where VFj_4(t) (resp. ^"4 (C/(r), A)) denotes the fourth component of lyj(T) (resp. T{U{t),X)). 
Hence, we have that: 

• Each solution U £ C°([0, ro],M'*) of (4.14) gives a solution of (4.7) that is bounded on [0,to]. 

• Every bounded solution U e C°{[Q,to],R*) of (4.7) satisfies (4.14) provided that we add d3V3(T) + 
d4V4{T) to the right hand side for an appropriate d £ M^. 

• Existence of solutions of (4.14) is given by the uniform contraction mapping principle for sufficiently 
small (^1,(^2) and A. 
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• The resulting solution U satisfies U{t) — J2j=i ^j^ji''') + <-^to(|^|M| + \d\^) on [0,to]. 

We also need to compute the quadratic coefficient in d in front of V4(to). Using a Taylor expansion, we find 
that this coefficient is given by 



1^ W4_4(s) (7^QQ(coshs)) ds ~ —V 
Jo 



' f°° 3 

/ (^ao(coshs)) sinhsds + o(l) 
Jo 



n 



4.3.2 The far-field equations and matching 

In this section, we look into the far-feild regime, where the radial variable t is large. We make spatial 
dynamical system (4.7) autonomous by augmenting the system with the equation 9^-6 ~ —e{2 + e) where 
e(r) = cothr — 1 to yield the new system 





i ^1 \ 


/ 


d 


U2 




d^ 


Us 






K^ ) 


\ 



-U2 



"3 

U4 

-Ml + M2 - (1 + f)u3 

(1 + e)?i4 + Alii + vui 
-e{2 + e) 



\ 



Wi 



(4.15) 



In the remainder of this section, we focus on the regime e 
denote A{oo, A) by the matrix 



^(00, A) 



which corresponds to the far field r ^ 1. We 



/o 





1 


0\ 











1 


-1 


1 


-1 





u 


-1 





-1/ 



where drU = A{t,X)U is the linearisation of (4.15) about the trivial state. We find that the matrix 
^(00,0) has four eigenvalues ao,aQ with multiplicity two (ag is defined in equation (4.8)); see Figure 15. 
As 3i{i') — —1/2, the trivial state f/ = is asymptotically stable at A = and then there is no bifurcation 
at the far field. Recall that in the Euclidean case a Turing instability occurs at infinity. In Figure 15, we 
summarize how the eigenvalues £ of ^(00, A) split close to A = 0. For A > 0, there exist four complex 
conjugate eigenvalues with di{£) = —1/2. For A < 0, there exist four complex conjugate eigenvalues with 
^{£) =^ —1/2 and stable manifold W^ is the union of the stable fast manifold (which we denote by W^ ) and 
the stable slow manifold (which we denote by W^) corresponding to the fast and slow decay to the trivial 
state. 
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; X 
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'-1 
1 2 
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c 
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',-1 
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'-J. 

; 2 

X 

I 



c 
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x=o 



A>0 



Figure 15: Schematic splitting of the eigenvalues £ of ^(00, A) for different values of A. At A = 0, the 
multiplicity is two. Eigenvalues in the red box correspond to the stable fast manifold W^ (A). 

First we argue that the centre-unstable manifold W™ and stable manifold W^ should intersect. We have 
that Vi{t) and V3{t) decay like e^'^/^, while V2{t) and V4(t) decay like re^'^/^ as t — > 00. Hence the tangent 
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space of the stable manifold at {u, A) = is spanned by (Vi(r), V2(t), V3(r), V4(t)). On the other hand, we 
showed in Lemma 4.1 that the tangent space of the core manifold is spanned by Vi{t) and V2(t). Then 
these tangent spaces would intersect along the two-dimensional subspace spanned by Vi{t) and V2(t). 



;^/ 



^;^ 



need to find an explicit description of W^ . To do this, we use successive, well chosen change of variables to 
put (4.15) into normal form. We first define the linear change of coordinates 



U^A 




\0j 



B 



( \ 

2ao + l 

1 
Vao(2ao + 1)7 



c.c, 



or equivalently, 



'A' 



iui 



■,Va 



{-\ui - \U2- 



-| {\U2 +U4) - i^U2 



and C/ = (mi,M2,U3, U4) 



In these coordinates, the linear part of (4.15) becomes at A = 0, 



drA = 



1 e 

-2-2+^ 


V3 
2 


e%/3 
6 


1 e 

-2-2+^ 


"V3 
2 


eV3 
6 



A 



B 



1 



i)-^a 



^u-a 



i\/3 



6 



drB 

drf- = -e(2 + e). 
Lemma 4.2. i^ia; < m < 00, then there exists a change of coordinates 



B, 



A' 
B, 



[i + r(6)] 



B , 



0{{\\\ + \A\ + \B\){\A\ + \B\)), 



so that (4.16) becomes 



(4.16) 



(4.17) 



drA = 


^ H-i)-^ 


drB = 


V 2 2; 3 


dre = 


= -e(2 + e). 



h.o.t.. 



XA + h.o.t. 



(4.18) 



T/ie coordinate change is polynomial in {A, B, e) and smooth in A and T(e) = 0(e) is linear and upper 
triangular for each e, while (p{r) satisfies 



dr(t>ir) 



V3 



Oi\X\ + \e\ + \A\'), 0(0) =0. 



Note that at (a. A) = (0,0), the trivial state {A,B) = (0,0) is hyperbolic such that the higher order terms 
in equation (4.18) are exponentially small for r ^ 1 and A small enough and can be neglected. We can also 
directly solve the linear part of equation (4.18), for A < 0, to obtain 



1 



A{r) \ ^ _^^ 

b{t) ) yiiS^ 



qie 



926 



(4.19) 



We want to find solutions that have a finite energy density with respect to the hyperbolic measure, i.e. 
functions that are in L^(M+, sinh(r)(iT). This restriction implies that we need to track the stable fast 
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manifold W^ (A) of equation (4.19) which corresponds to eigenvalues i of .4(oo, A) with real part less than 
— I as shown in Figure 15. Thus, for each fixed tq ^ 1 and for all sufficiently small A < 0, we can write the 
T = To-fiber of the stable fast manifold W^ (A) of equation (4.19) near U ~ as 



WI^(A) \r- 




-To/2 



-My^(i + o.„(iA|))( ^° 




(4.20) 



for ^ e C. 

We can now finish the proof of Theorem 4.1. To do this, we need to find nontrivial intersections of the 
stable fast manifold W^ (A) with the centre- unstable manifold yV™(A). To this end, we write the expansion 
(4.13) for each fixed tq ^ 1 in the {A, B) coordinates and afterwards in the coordinates {A,B). Using the 
expansions of the associated Legendre functions given in Table 2 we arrive at the expression 




-ro/2 



.(^.„+$„) / fd,{l + Oil))+Tod,[-^^ + Oil)) 



-d,{^^ 



0{l)) + ^[-v{I + o{l))dl 



+e 






(4.21) 



We can apply the transformation (4.17) to equation (4.21) and obtain the expansion 

^A'^0+O[r-^)+O^„{\\d\ + \d\^) ( 0^0(^1^1 + 1^12) \ 

l^ OrMd\ + \d2? + \di?) ) 

^di(l + 0(1)) + rod2(-i^ + 0(1)) 



>V™(A) 




_^ gi(<['0+O(r-^)+O,„(A|<J| + |d|2) 



-d2(*^+0(l)) 

\ 



^(-^(I + o(l))J?) 



(4.22) 



The final step of the analysis consists in finding nontrivial intersections of the stable fast manifold W^ (A) 
given above in equation (4.20) and the core manifold yV™(A) given in (4.22). We can easily solve this 
problem in (^1,^2) to find that 

2x/^ 



rfi = 



vCiI 



, and 6^2 = 0(A). 



Then c = 



vCiI 



in equation (4.6) and sign(c) — sign(zy). This completes the proof of Theorem 4.1. 



4.4 Numerical computation of spots 

In this section, we describe the use of numerical continuation (and the continuation package AUTO07P [26]) 
to compute solutions of the systems of ODEs described by (4.7). Solutions of these spatial dynamical system 
correspond to steady states of the Swift-Hohenberg equation (1.1) where the radial coordinate r has been 
recast as time in AUTO07p's boundary value problem (BVP) solver. The BVP is set up on the domain 
T g [0, L] with homogeneous Neumann boundary conditions given by ^2(0) = 1*4(0) = ui{L) = U3{L) = 0. 
Typical parameters for the AUTO07P radial computations are L = 1000 and AUTO07p's NTST=400 with 
standard relative tolerances that are specified in AUTO07p's manual. Initial data for continuations are 
obtained by first solving the radial Swift-Hohenberg equation (4.4) with the fsolve routine of Matlab and 
then using a parameter continuation to add in the radial terms. 

In Theorem 4.1, we have shown the existence of spots of equation (4.4) for any fixed v j^ and ry e M when 
< |A| <C 1. We compute spot solutions for i^ — 0.8 and rj — 1 and summarize the results in Figure 16. Spots 
do bifurcate off u = at A = and turn around at a saddle-node bifurcation at A = —0.0425. At this fold, 
spots regain stability with respect to radial perturbations, but they remain unstable with respect to general 
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Figure 16: Spots of (4.4) are shown for v = 0.8 and rj = 1. (a) Bifurcation diagram of spots, (b) Log-plot 
of the maximum height w(0) versus the parameter A: the sohd Hne is the result of numerical computations, 
which we fit with the dashed line 0.531og(— A) + 0.87. (c) Profile of a typical spot, (d) Decay rate of the 
tails of the solutions rescaled solutions M(r)e'^/^ in black at A = —0.015. Linear best fits is shown in gray. 

perturbations. The corresponding bifurcation diagram (see Figure 16(a)) is similar to the one computed 
in Euchdean geometry by Lloyd & Sandstede [51]. The computations confirm the scaling u(0) w VlAJ as 
A — > 0. We also verify in Figure 16(d) that the spots found in Theorem 4.1 are L (M+,sinh(T)(iT) functions. 
Close to onset at A = —0.015, we plot the decay rate of the tail of the rescaled spots u(T)e'^/^ which confirms 
that spots decay faster than e^"^/^ and hence are in L^(M+, sinh(r)(ir). 

Now, we compute spots for v — 1.6 and -q — 1 and show the results in Figure 17. Panel (a) shows a bifurcation 
diagram in A where the branches are represented in terms of the Euclidean L -norm |Jm|| = (/p°° uItYtcIt) 
This is due to numerical difficulties of computing the hyperbolic L norm presented by the sinh function in 
the integrand but, as indicated in [33], the Euclidean radial norm is a good solution measure. The panels 
(b)-(e) show solution profiles at different points on the bifurcation diagram for fixed value of A as indicated. 
As one moves up on the branch, rolls are added one by one to the tails of the spots: this corresponds to 
adding concentric rings that surround the spot. Panels (b)-(e) show that the amplitude at the core is still 
much larger than the amplitude of the concentric rings that are added. The bifurcation diagram in panel (a) 
presents similar characteristics of "snaking" -type diagram found in Euclidean geometry with what could be 
an analog of a Maxwell point at \m — —0.2218 [6, 11, 13, 51]. Note that this "hyperbolic snaking" structure 
was not reported in [33] in the case of neural field equation. 



5 Horocyclic traveling waves 

The basic material of this section was introduced in [17] in another context (analysis of a model equation 
for the detection of textures by the visual cortex of mammals) . 

Instead of being periodic on a lattice in D, we may look for states which assume the form of hyperbolic 
plane waves (or horocyclic waves) as defined in 3.1.1. Let us consider the horocycles which come in contact 
with dn at the point h — \ (see Figure 3). Horocyclic waves with base point b would be constant along 
the horocycles and periodic along the "hyperbolic" coordinate. In other words, writing z ~ Ug ■ Ot ■ O in 
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Figure 17: Bifurcation diagram in A of spots of (4.4) is shown in panel (a) for v = 1.6 and -r] — \. Panels 
(b)-(e) contain plots of spots at A = —0.2218 at different values of ||w||. 

horocyclic coordinates, the solution would satisfy the invariance properties 

uijis ■ Ut ■ O) — ui^a-r ■ O) and u[a.r-\-T ■ O) — u{ar ■ O) 

for some period T. Such solutions would be the counterpart of the "stripe" or "roll waves" solutions which 
occur for the Swift- Hohenberg equation and most equations which allow for pattern selection in the Euclidian 
plane. We shall see in this section that such solutions do indeed bifurcate in D and have many common 
features with the Euclidean case with some important differences. In particular, these waves are traveling at 
constant, generically non zero speed in D, while they are steady in the Euclidean plane. 

Assuming horocyclic invariance as defined above reduces the coordinates to the single variable r. The 
Laplace-Beltrami operator then reduces to (see [38] where the normalization of the transformations ar is 
slightly different from ours) 

A" 






d_ 



Now the equation (1.1) reads 



Ut 



-{a^ + lS!^fu + \u + vu^ - Tju^, 



z e 



(5.1) 



and w is a function of r and t only. 

As was mentionned in 3.1.1, the linear stability analysis of the trivial state of (1.1) with respect to the 
elementary eigenf unctions Cp ;, of the Laplace-Beltrami operator comes back to solving the "dispersion re- 
lation" (3.1). These "hyperbolic plane waves" are precisely functions of the coordinate t only and they 
are in exact correspondance with the elementary eigenfunctions of Ap, which have the form e*'^'^, A; e C 
The correspondance is given by the relation ik — ip + 1/2, which translate for the eigenvalues of Ae to 
— fc^ — ik — —p^ — 1/4. In order for these "waves" to be periodic in r we require fc e M. Observe that the 
eigenvalues are still complex. The dispersion relation (3.1) for perturbations e'^*+*'='^ now reads 

o = -k:^ + (2a2 -I- l)k^ - a'' + A + 2k{a^ - k^)i 

With no loss of generality we now assume that a^ = 1/2. Then from the dispersion relation it follows that 
the most unstable modes occur at k — ±1 with critical parameter value A = —3/4, which corresponds to 
eigenvalues a — ^i. Therefore the bifurcation of periodic hyperbolic plane waves is a Hopf bifurcation to 
time periodic solutions. 
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This is however a simple case of Hopf bifurcation thanks to the translational invariance of Equation (5.1). 
Indeed let us look for solutions in the following form: u(r, t) = v(r— (l+w)i) = v{y) with lo to be determined. 
This is a uniformly traveling wave at constant speed 1+lu. Then (5.1) becomes 







1 d^ 



dy' 



fl 



- (jj) \- AV + vv 

dy 



rjv 



(5.2) 



At A = —3/4 and a; = the linear part of this equation has an eigenvalue at with eigenvectors ^o = e*** 
and Co and solving this bifurcation problem is a classical exercise. 

We set A = —3/4 + /x and v[y) = AC^q + AQq + '^{A, A, w, fi) where ^ is a map on the orthogonal complement 
of the kernel of the critical linear part spanned by Coi Co (we consider here the usual L^ inner product of 
27r periodic functions). The Taylor expansion of ^ and the bifurcation equation for A, A and u) are then 
calculated using the Lyapunov-Schmidt decomposition method. 

The following equations are obtained at leading order: 



= (^ - iuj)A + 2v 



78 -4i 
75 



-■ir]]A^A + A-o{\\A\\ 



m 



and complex conjugate. The bifurcated solutions are therefore given at leading order by the amplitudes 



1^11' 



156 
"75" 



v^ -2,'n 



M 



and frequency 1 — 8/75t^^||74|p. Observe that the bifurcation is supercritical if 77 > and v is small enough, 
and that the speed of the traveling wave is lower that the critical one ii v ^ Q (it depends on higher order 
terms if i^ = 0) . 





(b) 



Figure 18: Left: Color plot of a solution of (5.1) in horocyclic coordinate r. Right: Snapshot of the same 
solution profile of (5.1) in the Poincare disk at time t = 200. Values of the parameters are: /i = —0.01, 
V = 2.7 and r/ = 1. 

In Figure 18, we show the example of a bifurcating horocyclic traveling wave of equation (5.1). Numerical 
simulation of traveling wave is carried out using MATLAB. We take as initial condition a small localized 
solution around r = and run the simulation for t G [0,200]. We use a semi-implicit finite differences 
method to compute the solution of (5.1). Space and time discretizations are taken to be At — 0.05 and 
At = 0.02. The values of the parameters are taken to be /i = —0.01, v = 2.7 and ?] = 1. It can be seen 
from Figure 18(a) that the solution converges to a traveling wave solution in horocyclic coordinate. We also 
plot in Figure 18(b) the corresponding solution profile at time t = 200 in the Poincare disk where the white 
arrow indicates the direction of propagation of the wave. 
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6 Discussion and some open problems 

In this review, we have analyzed the bifurcation of i) spatiahy periodic solutions, ii) radially localized solutions 
and iii) traveling waves for the Swift-Hohenberg equation with quadratic-cubic nonlinearity defined on the 
two-dimensional hyperbolic space: D (Poincare disk). 

For the bifurcation problem of periodic solutions, we have made use of the concept of periodic lattice in 
D to further reduce the problem to one on a compact Riemann surface D/F , where F is a co-compact, 
torsion-free Fuchsian group. Following the method introduced by [18] in a different context, we have applied 
techniques from equivariant bifurcation theory in the case of an octagonal periodic pattern, where we have 
been able to determine the generic bifurcation diagrams of each irreducible representation. As a case study, 
for the irreducible representation associated to the lowest nonzero eigenvalue, we completely computed the 
coefficients of the reduced equation on the center manifold and determined in parameter space {v, fi) the 
stability of each bifurcated branch of solutions. 

To prove the existence of a branch of bump solutions bifurcating near onset we have followed the presentation 
in Faye et al. [33], which applied techniques developed by Scheel [60] and Lloyd & Sandstede [51] . In the 
Poincare disk, the analysis near the core manifold required the development of a detailed knowledge of the 
asymptotics of the associated Legendre functions, which are the counterpart of the Bessel functions in the 
Euclidean case. It turns out that the essential difference between the two geometries comes from the far field. 
At infinity, Bessel function Jo(?') scales in term of the radial coordinate r as l/y^ whereas the associated 
Legendre function V_x .^{cosh.{T)) scales in term of the polar geodesic coordinate r as e^'^/^ for a = 1 
in equation (1.1). Moreover, in the Euclidean case and for the trivial state at infinity, there is a Turing 
instability bifurcation. However, in the hyperbolic case, the trivial state is always asymptotically stable at 
infinity and this simplifies the resulting matching problem. We have also used numerical continuation to 
track branches of radially symmetric solutions away from onset and found for high value of the parameter v 
what could be the analog of a snaking diagram for the Poincare disk (see Figure 17(a)). 

Finally, we have seen that a branch of traveling waves invariant along horocycles and periodic in the transverse 
direction, bifurcates off the trivial state through of Hopf bifurcation. This striking result highlights once 
again, the inherent differences that exist between pattern formation in Euclidean and hyperbolic geometry. 

This analysis is far from being complete and there still exist many open questions. We list below some 
problems that we consider to be important. 

• We studied the bifurcation of H-planforms for the regular octagonal lattice. There are infinitely many 
types of periodic lattices in the Poincare disc. Of course, the methods presented in Section 3 would 
easily be adapted to other lattices. However this raises the question of the observability of such patterns 
in a natural system or under direct simulation of the Swift-Hohenberg equation. Indeed, in addition 
to the high degeneracy of the bifurcation problem, there is a rigidity of the lattices in hyperbolic space 
which, as explained in Section 3, could make the observation of such patterns unlikely. We may imagine 
mechanisms such as spatial frequency locking which could overcome this difficulty, but this is a very 
challenging problem. 

• In section 4, we presented an existence result of localized radially symmetric solution. Other localized 
radial states can be expected and the theory developed in [51, 53, 55] could be extended to yield the 
existence of rings and other types of localized radially symmetric solution, which were not investigated 
in this paper. We also think that the existence of horocyclic traveling waves might lead to the formation 
of target or source emitters. 

• In the Euclidean case, localized hexagonal patches have been studied in the Swift-Hohenberg equation 
with quadratic-cubic nonlinearity [52] . It would be of interest to numerically investigate if one can find 
localized patches with Dg-symmetry for the range of parameters (y, /i) in Figure 14 where the branch 
of solution with Dg-symmetry is stable. Note that Faye et al. [33] have numerically computed branches 
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of solutions with Dg-symnictry bifurcating off the branch of radially localized solutions for a neural 
field equation set on the Poincare disk. 

Acknowledgment: GF is grateful to James Rankin and David Lloyd for their helpful comments on the 
functionahty of AUTO. 
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